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ABSTRACT 



This thesis demonstrates a method to determine optimum, tactical movement paths 
for a specified vehicle and or small unit based on the operator's cognitive decision 
processes, as well as the physical effects of terrain and environment on mobility. The 
approach uses psychometric techniques inherent to the Generalized Value System (GVS) 
to determine a “Power Function” based on a specific tactical scenario and given equip- 
ment configuration and provides a means to determine the Tactical Movement Potential 
(TMP) for each terrain cell. This cognitive value in an interval scale can then be trans- 
lated into the same scale as the physical continuum using techniques proposed by L.L. 
Thurstone and W.S. Torgerson. The cognitive time value based on the user's decision 
process is then added to the physical traversal times for each cell computed from output 
provided by The Condensed Army Mobility Management System (CAMMS). This ren- 
ders a value mapping which can be optimized using one of several existing algorithms. 
The Dijkstra Algorithm is used in this demonstration model. The resulting sets of path 
points are optimized for speed time and the cognitive tactical considerations evaluated 
using these psychometric methods. T he movement path and resulting times can be used 
in combat planning and modelling. This output is also particularly important in deter- 
mining the time values needed to compute the Situationally Inherent Power (SIP) of the 
GVS. This methodology could be applied to almost any tactical decision process in the 
development of expert systems and models. 
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THESIS DISCLAIMER 

The views expressed in this thesis are those of the author and do not reflect the of- 
ficial policy or position of the Department of Defense or the L'.S. Government. 

The reader is further cautioned that certain vehicle system input parameters used 
and portions of the computer program developed in this research are not valid for all 
scenarios of interest. While every effort has been made, within the time available, to 
verify all computer programs: they cannot be considered validated. Any application of 
the methodology using these programs is done so at the user's own risk. 
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I. INTRODUCTION 



A. GENERAL 

The ability to accurately portray the battlefield environment is a key requirement in 
combat model development. This same ability is required in the construction of tactical 
decision aids and tools. One major factor to be considered in any of these applications 
is the effect of terrain, vegetation, and weather on the mobility trafiicability of weapons 
systems and units. Another major factor is the effect of these same parameters on visi- 
bility, detection, and acquisition by various weapons systems. While the preceding fac- 
tors are objective for the most part and can be accurately measured; other factors are 
not so clearly defined and may possess large variances. This is certainly the case where 
a human decision process is required such as the decision logic used by various combat 
unit and vehicle commanders in tactical movement. A clear distinction must be made 
between mobility as modelled in relation to the physical constraints of the battlefield and 
movement which requires interaction with the user decision maker or must be modelled 
taking the pertinent cognitive considerations as well as the physical constraints of the 
battlefield into account. 

In any case, all of these factors and more must be accounted for in our combat 
models and tactical decision tools. There are a number of existing high resolution 
combat models and tactical decision aids which incorporate many of these functions. 
Most of these packages use one of several available versions of digitized terrain data and 
either actual, historical, or user input weather data. Most combat models also incorpo- 
rate appropriate tactical decision logic algorithms to control some key functions such 
as detection, acquisition, and other essential combat actions of player units. Some tac- 
tical decision tools also do this. However, at the present time there is no combat model 
or tactical decision tool which combines the appropriate decision logic to determine the 
optimum path a unit or vehicle should traverse in order to simultaneously maximize 
speed and trafiicability, minimize inter-visibility and detection, and appropriately bal- 
ance their effects in terms of the decision logic. 

The need and desire of commanders and combat planners to have and use such 
tactical decision tools and to incorporate them into combat models is of definite benefit. 
Some commanders presently have a tactical decision aid in the form of the Army Mo- 
bility Model (A.MM) and the Condensed Army Mobility Management System 
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(CAM MS) which provide trafTicability. speed, and time contour mappings for use in 
route selection and other tactical decisions. There are also any number of combat 
models available for study where a commander or combat analyst can select appropriate 
terrain data bases and manually input desired unit or vehicle path points and or posi- 
tions. This is the case with present JANUS based models, the Battalion Combat Out- 
come Model (BCOM), and others. However, there is no tactical decision tool or combat 
model that can systemically determine sets of optimum, tactical path points for the 
commander or for use in analysis. This is also true of almost all models where tactical 
decision processes are represented. 

The first requirement in developing such a model is to properly identify those vari- 
ables which could and or would affect the tactical movement or other decision process 
such as the mission, threat, corresponding locations, ranges to various threats, vehicular 
agility or localized speed, allowable completion time, and others. It is soon obvious that 
a large number of variables, many having probable interactions with other variables, 
exist. The possible combinations or instances is quite large and each variable may vary 
significantly between instances. Structuring this maze in order to be able to extract usa- 
ble data on the cognitive process is no small task. 

Another requirement is finding a theoretically sound method for evaluating large 
groups of variables simultaneously and then being able to compute a specific value 
which appropriately measures their combined effect on the process. This field of study 
has a name all its own. Psychophysics. 

A third requirement is to find a theoretically sound method that allows for the 
mapping of these psychophysical values into an appropriate physical continuum so that 
they may be used concurrently with physical attributes and variables in the model. This 
field of study also has a name all its own, Psychometry. 

The art of combining and applying these techniques is known as Psychometrics. 
We will use these techniques in conjunction with an existing mobility model, CAMMS 
(Version 2.0), to provide total movement time values based on physical engineering 
traversal times and the translated psychometric times for each terrain cell and subse- 
quent traverse segment on the battlefield. These values will then be used as network cost 
values in a form of Dijkstra's Algorithm which renders the set of minimum cost optimum 
path points for that scenario. 

The organization, purpose, scope, and assumptions of this research are explained in 
this Introduction. Background information on previous route selection methodologies, 
terrain representations, physical attributes of terrain considered in mobility, cognitive 
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considerations which impact on movement, and the desire to place such technologies in 
the hands of the commander in the field is presented in Chapter II. The development 
of the “Psychometric Function” using psychometric procedures inherent to the Gener- 
alized Value System, GYS, is shown in Chapter III. Application of these techniques and 
construction of the model algorithm is detailed in Chapter IV. Results of the algorithm 
implementation on an AT-3S6 personal computer and limited verification of these results 
are presented in Chapter V. Proposed enhancements to the existing model and areas re- 
quiring additional study are presented in Chapter VI. Finally, Chapter VII discusses 
possible utilizations for this methodology and proposed applications with respect to 
tactical decision and training tools, research and analysis tools, combat models, and the 
Generalized Value System. 

It is of significant interest to note that prominent Soviet military analysts. Nikita 
Moiseyev in particular, have proposed a model framework based on decision theoretic 
developments using such state variables and decision functions since the 1970's. Two 
specific areas discussed have been tactical decision processes and movement of the line 
of contact. [Ref. 1: pp. 21-26] 

B. PROBLEM STATEMENT 

1. Purpose 

This thesis examines and demonstrates a methodology for determining opti- 
mum. tactical movement paths for a specified vehicle and or small unit based on the 
operator's cognitive decision processes, as well as the physical effects of terrain and en- 
vironment on mobility. The resulting sets of path points are optimized for speed time 
and the cognitive tactical considerations evaluated using psychometric methods thereby 
rendering a truly optimum, tactical movement path and optimum traversal time. The 
movement path and resulting times are to be used in combat planning and modelling. 
This output is also particularly important in determing the time values needed to com- 
pute the Situationally Inherent Power (SIP) of the GVS. 

2. Scope 

This research will be limited to one specific scenario of interest in order to 
demonstrate the methodology without undue involvement in various scenarios with dif- 
ferent sets of fixed state variables. The specific scenario used is discussed in detail in 
Chapter III. This research will not involve detailed studies of relative effective range used 
in the algorithm and discussed in Chapter IV. There will also not be any extensive ver- 
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ification and validation of the results of the computer application of the algorithm due 
to a lack of time and resources. 

This research will examine the use of a movement model based on the method- 
ology for a single vehicle unit of specified type and mission against a threat array of 
specified types performing a certain mission while in a selected area of operations and 
known terrain. 

3. Assumptions 

1. The mathematical forms of range dependency equations studied by Bonder and 
Farrell and used to portray perceived lethality as relative effective range curves re- 
lated to those found in JANUS are correct. 

2. The survey results, interval scales, and regression analysis are valid for the range 
of values surveyed only. However, these should encompass the vast majority of 
instances. 

3. The translated cognitive time scale and traversal time scales are additive and not 
multiplicative. 

4. This methodology is applicable to all types of tactical decision processes when the 
appropriate state variables can be identified. 
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II. BACKGROUND 



A. PREVIOUS ROUTE SELECTION AND TACTICAL MOVEMENT MODELS 

1. The Dynamic Tactical Simulation (DYNTACS) Model 

The DYNTACS Model was developed by the Systems Research Group at Ohio 
State University from 1964 to 1969. The model is extremely high resolution in order to 
predict individual unit performance in armored combat engagements from both the de- 
sign and operational perspectives (Ref. 2]. Dynamic programming is used to run and 
solve the simulation which consists of a driver program and 34 subroutines. At least 15 
of these algorithms are related to or affect tactical movement in the model [Ref. 3]. An 
explicit surface terrain representation of the battlefield is used in order to determine 
plane departure points for a predetermined grid size [Ref. 3: pp. 57-66]. This represen- 
tation is closely related to and can be derived from explicit grid terrain data available 
from a number of sources [Ref. 4: pp. 3.3-3.15]. 

The model determines an optimum route for an advancing element or unit by 
optimization where “Tactical Difficulty" is minimized. Tactical difficulty is computed 
using a heuristic of the form, TD = (1 + E)T , where E is the difficulty computed due to 
a set of factors along each route segment and T is the travel time for that segment [Ref. 
3: pp. 94-116]. Such heuristic procedures are computationally efficient but do not guar- 
antee an optimal solution. The set of factors and corresponding values for the difficulty, 
E. were obtained using comparative judgement techniques. Travel times are computed 
using engineering models for each weapon system. It is interesting to note that difficulty 
is a function of the effective range, actual range, and threat disposition [Ref. 3: pp. 
97-100] and that psychometric methods similar to many we will examine were used to 
obtain these values. The major shortcoming with the DYNTACS approach is that the 
cognitive and physical scales are related multiplicativelv using a heuristic due to existing 
computer limitations rather than developing a truly optimal technique. 

2. The Simulation of Tactical Alternative Responses (STAR) Model 

The STAR Model is similar to DYNTACS in most respects except that dynamic 
programming is discarded and the optimization is performed in FORTRAN applying 
Dijkstra's Algorithm. The additional concept of the “Sliding Pattern" is used to restrict 
the search for the path [Ref. 5: pp. 23-34], The route selection is performed sequentially 
from start point, to a horizon or horizons, to an objective by sliding the optimization 






along the network. Unfortunately, the identical heuristic equation found in DYNTACS 
is also applied in this model [Ref. 5: p. 33] rendering it suspect. 

3. Other Route Selection Methodologies 

There are any number of similar movement models proposed over the past 25 
years and certainly all of them cannot be mentioned. However, it is necessary to intro- 
duce another type of model which is probabilistic in nature. One such model was inves- 
tigated and formulated by Sitmourang [Ref. 6] in 1981. This dynamic route selection 
model performs an optimization in order to minimize the probability of being attrited 
enroute to the destination. The attrition probabilities along each route segment are de- 
rived from two sources. These are enemy elements capable of destroying the tactical el- 
ement and terrain factors which end movement. Speed reduces the probability of 
attrition due to the enemy but increases the probability of attrition due to terrain, espe- 
cially rough terrain [Ref. 6: pp. 21-22]. The total probability ofbeing attrited. P. is taken 
to be the sum of the probability ofbeing killed by the enemy, PKFOE, and the proba- 
bility ofbeing killed due to terrain speed, PKVEL, which is minimized using Dijkstra's 
Algorithm [Ref. 6: p. 34]. 

There are a number of problems with this methodology. The model assumes 
perfect information in computing probability of kill due to enemy. This is certainly not 
the case in most tactical decisions. The probability of attrition due to terrain speed is 
computed using a heuristic equation with a power factor or constant for all combina- 
tions of conditions which is applied to a ratio of recommended speed to actual speed. 
This value is then used multiplicatively with normal expectancy in order to determine a 
probability of attrition due to speed in given terrain [Ref. 6: pp. 24-26]. No reference or 
theoretical support for this technique is provided. This certainly makes such a method- 
ology suspect. 

B. CURRENT MOBILITY MODELS 

Mobility models provide comprehensive mappings by terrain cell or other construct 
of the measure or ability of weapons systems to physically traverse types of terrain in 
varying conditions. Mobility models should not be confused with movement models 
which attempt to describe how tactical elements traverse the battlefield based on both 
physical and cognitive considerations. The two primary mobility models used today are 
The NATO Research Mobility Model (NRMM) and The Condensed Army Mobility 
Management System (CAMMS). Most mobility routines used in models today are based 
on this family of models. The NRMM and CAMMS are both derived from the Army 



6 



Mobility Model (AMM) developed in the 1970's. The CAM MS is currently used as the 
base model for: 

1. The Planning Analysis Work Station (PAWS) 

2. The Terrain Analysis Work Station (TAWS) 

3. The Airland Battlefield Environment (ALBE) 

4. The Digital Terrain Support System (DTSS) 

The model is also used by numerous Army, Defense, and other government activities to 
include The Defense Mapping Agency, The L’.S. Army Command and General Staff 
College, The United States Military Academy, The Army Research Institute, The U.S. 
Army Engineer School. The U.S. Army Training and Doctrine Command (TRADOC). 
and The TRADOC Analysis Command (TRAC). The CAMMS is used as the base mo- 
bility model in this research due to its applicability and acceptance in the modelling 
community. [Ref. 7] 

1. The Condensed Army Mobility Management System (CAMMS) 

The CAMMS vehicle performance model was originally developed by simplify- 
ing the existing AMM through addition of a vehicle preprocessor and restricting the 
model to one vehicle type at a time. Obstacle performance is also determined statistically 
based on years of testing and data collection. This allows the CAMMS to operate on 
personal computers versus larger mainframe computers for which the AMM was devel- 
oped. The model is much too large and detailed to be discussed fully here; however, an 
examination of the physical aspects considered in the CAMMS, output available, and 
current uses of the model is worthwhile. 

a. Physical Considerations Affecting Mobility 

The CAMMS off road prediction model evaluates three general categories 
of terrain data. These are surface geometry data, surface composition data, and vege- 
tation data. Surface geometry data includes slope information, surface roughness factors, 
and all obstacle data. Surface composition data primarily relates to the type of soil and 
its strength which are greatly affected by weather. The vegetation data includes stem 
diameter, type, and spacing data along with visibility information. These data are used 
in a number of routines to compute surface traction, resistance, obstacle traction, inter- 
ference, gap crossing capability, vegetation effect, and ride dynamics which are combined 
to produce a highly accurate speed prediction for the vehicle. [Ref. 7] 
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b. CAMMS Output Mappings 

The CAMMS is widely used throughout the Army and particularly in the 
European Theater of Operations. A range of mappings or Tactical Decision Aids. 
TDA's, are currently available in the CAMMS. These include: 

1. Vehicle Predictions 

a. Speed 

b. Time 

c. Reason Limited 

d. Comparison of Two Vehicles 

2. Foot Soldier Prediction 

3. Maneuver Damage Prediction 

4. Soil Strength Prediction 

5. Route Analysis Prediction 

The vehicle predictions and especially speed and time mappings for the digitized terrain 
grid are of most significance to this research [Ref. 7). These will be addressed concerning 
their potential use in subsequent chapters. 

c. Technology to the Field 

The CAMMS takes advantage of recent advancements in the personal 
computer field to provide an extremely powerful yet portable and inexpensive tool for 
the commander in the field. The most current CAMMS requires only an IBM compat- 
ible AT personal computer with a math coprocessor. DOS 3.2, EGA graphics, 1 mega- 
byte of RAM, and sufficient hard disk drive storage in order to operate. All the TDA's 
mentioned previously are available in a matter of minutes provided the appropriate 
digitized terrain data bases are available. The only shortcoming in the current system is 
the lack of standarization in terrain data bases and lack of data in general for some re- 
gions. The CAMMS is being used extensively as an aid on REFORGER exercises and 
as an operational planning aid in some tactical units. 

The methodology developed during this research will also be derived so as 
to make it conformable for use on a personal computer. This will allow the model to be 
used as a portable tactical decision tool in the hands of the commander. 

C. COGNITIVE CONSIDERATIONS IMPACTING ON MOVEMENT 

The discussion up to this point makes it quite apparent that tactical movement 
modelling is much more involved than mobility modelling. The major question becomes, 
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how do we model the cognitive decision process? There is no definitive answer to this 
question and numerous techniques are available and feasible depending on the applica- 
tion. It is apparent that the mission, terrain effects, and threat all influence the 
combatant or tactical decision maker to some degree. The decision concerning where and 
when to move is made based on the combatant's evaluation of how the enemy can affect 
his mission accomplishment and or survival as he moves to each possible location. There 
is also a question of the trade-off made by most humans between surviving and accom- 
plishing a mission within a certain time. Some people are more risk averse than others. 

The process is synonomous to shopping for fruit or vegetables of different types and 
deciding which items to purchase based on their size, shape, color, texture, and price. 
The shopper makes a determination when comparing different items as to which is best 
in terms of overall quality relative to all others and what he will purchase with a certain 
amount of funds just as a combatant makes a determination of the potential lethality 
facing him based on the types of enemy, their range, and their capabilities based on 
range and disposition in order to decide what risk he can accept when moving on the 
battlefield. 

An enumeration of the variables identified as relative to this tactical movement de- 
cision process and how they are structured for use in this research is provided in detail 
in Chapter 111. 
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III. PSYCHOMETRICS AND PSYCHOPHYSICAL ANALYSIS IN THE 

CVS 



A. GENERAL 

The Generalized Value System, GVS, is an axiomatic value system used as a tool for 
evaluating the power or potential of entities at present and future times. The formal 
structure of GVS was first proposed by Professor Arthur Schoenstadt of the Naval 
Postgraduate School in 1985 [Ref. 8: p. 3). This procedure outlines a method of deter- 
mining the state of an entity over time, allows us to compare the values of entities rela- 
tive to one another, and is consistent for all entities on the battlefield. The values 
assigned to each entity are referred to as “power” in previous GVS documentation. 

Power is a subjective value given to each element on the battlefield and combines 
both inherent and derived power. Derived power consists of the power an element pos- 
sesses due to its ability to influence support other friendly elements [Ref. 9: p. 33]. In- 
herent power is the ability of an element to directly affect the outcome of a battle [Ref. 
9: p. 32]. Inherent power may be categorized as Basic Inherent Power (BIP) or 
Situational Inherent Power (SIP) which are the power an element has at full strength 
and the power an element is predicted to have at some future time, respectively. The 
SIP is highly dependent on the available time. t u , which is determined by the time do- 
main networks for the respective battlefield [Ref. 8: p. 5]. 

The methodology proposed in this research directly supports use of the GVS by 
providing a means of determining actual optimum time values for the terrain grid net- 
work representing the battlefield. The methodology will also provide actual optimum 
time values and sets of path points for use in other combat models using conventional 
Lanchester attrition or firepower scores as long as they utilize explicit terrain grid net- 
works for representing the battlefield. It will also provide excellent tactical movement 
planning and decision tools when incorporated into the appropriate mobility models. 

Modelling the cognitive process in order to determine appropriate values for a gis'en 
tactical movement scenario and instance and to then be able to map these values into 
the physical continuum, time, on the battlefield requires an understanding of how atti- 
tudes are measured, how to develop respective interval scales, how to convert cognitive 
values in interval scale into physical values in a ratio scale, and how to solve for the 
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“Psychometric Power Function" which allows the user to solve for this value given any 
set of realistic variables at different levels. 

B. THE MEASUREMENT OF VALUES 

The correct measurement of psychological values is fundamental to evaluating in- 
stances or combinations of variables and subsequently being able to use them in any 
physical construct. However, there are numerous ways to develop measurements and 
each results in different scales with varying characteristics. The three major character- 
istics of measurement are order, distance, and origin. Distinctions between types of 
measurement and types of scale depends on what combinations of these characteristics 
exist [Ref. 10: p. 15]. 

A nominal scale exhibits none of these characteristics and will therefore not be 
evaluated in terms of potential use [Ref. 10: p. 17]. An ordinal scale does possess order 
and may have an origin; however, the scale does not exhibit distance which results in 
no ability to discriminate between instances in terms of a value [Ref. 10: p. 16]. This 
limits potential types of usable scales in psychophysical analysis to interval scales. 

1. Interval Scales 

Interval scales always possess the properties of order and distance. It is also 
often possible to determine a rational origin for such scales [Ref. 11: p. 196]. Interval 
scales which exhibit all three major characteristics are known as ratio scales. These scales 
are the most important scales used in math and science and the term “measurement” is 
often restricted to these types of scales [Ref. 10: p. 31]. Interval scales can be determined 
in a variety of ways for sets of instances and scenarios. Three powerful methods are 
Continuous Response Scale Judgements, Comparative Judgements, and Categorical 
Judgements. 

a. Continuous Response Seale Judgements 

A continuous response scale allows a judge to indicate the value of an in- 
stance on a continuous scale between a minimum and maximum bound, normally 0 to 
100. This method lends itself to Analysis of Variance (ANOVA) testing commonly used 
to determine relationships. The scale value of each instance is easily obtained and the 
summary statistics along with the distribution of responses for each instance are also 
easily obtained. However, one problem with this approach is that judges have difficulty 
actually determining values to assign between various instances [Ref. 8: p. 10]. This is 
magnified the larger the number of instances becomes. For example, only live variables 
broken down into five distinct levels would result in 5 s or 3125 unique combinations of 
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variables requiring responses. It is extremely difficult for a judge to realistically place a 
value on each with this number of instances and this approach was rejected due to this 
factor. 

b. Comparative Judgements 

The Law of Comparative Judgements presents the judge with all possible 
pairs of instances in order to determine which has the greater value in each pair. This 
discriminal process relies on the principle of “just noticeable differences" developed by 
L.L. Thurstone. The results of these judgements can be tabulated to provide a discrimi- 
nal difference and deviation [Ref. 11: pp. 39-47]. Application of this theory is done using 
the method of paired comparisons. A matrix solution based on the frequencies of the 

responses is used to compute an interval scale. However, with n instances it requires 
n{n — l ) 

- comparisons [Ref. 10: pp. 166-179). For the five variables broken down into 

five levels, this results in 3125 instances and 4,881.250 separate comparisons. This is far 
more comparisons than is feasible in most studies and too labor intensive in terms of 
data collection. 

c. Categorical Judgements 

The Law of Categorical Judgements parallels the Law of Comparative 
Judgements in many respects. The judges evaluate each instance independently and se- 
lect a rating category they feel best represents the general value of that instance. The 
categories are understood to be a mutually exclusive set of successive intervals which 
collectively exhaust all possible values of the continuum. Each rating category has de- 
scriptors and explanations which serve to help the judge with his task [Ref. 12: p. 1]. 
Normally distributed responses for each instance across the sampled population is a key 
assumption of this methoJ. This is easily observed by examining the frequency matrix 
in Appendix B. Computing the interval and or ratio scale is relatively straight forward 
using a matrix solution [Ref. 10: pp. 205-239]. A clear, concise method is given by Pro- 
fessor Glenn Lindsay. Discussions with Drs. Parry and Lindsay at The L'.S. Naval 
Postgraduate School assured that such an approach would be theoretically sound and 
yield usable results without requiring excessive survey and data collection effort and time 
once the scenarios were properly structured. 

C. THE EXPERIMENTAL DESIGN 

As previously addressed in the Problem Statement, the purpose of this study is to 
develop and demonstrate a methodology for determining an optimum, tactical move- 
ment path based on both the cognitive and physical considerations or constraints of the 
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battlefield. As was also stated in the Introduction, devising and structuring tactical sce- 
narios which insure the correct combinations of state variables or instances impacting 
on movement logic are addressed to the appropriate sample population can be an enor- 
mous and complex task. Administering the data collection effort once scenarios are de- 
veloped is another major step in the experiment. 

The process becomes more structured once the data are obtained. The Law of 
Categorical Judgements is used in this case and a well defined method for obtaining and 
translating the interval scale from the responses is available. Multiple regression analysis 
can then be applied to determine the functional relationship of the scale values for each 
instance to the set of state variables. The resulting regression equation becomes the 
“Psychometric Power Function” needed to determine values for use in GYS and other 
models. 

1. Scenario Development 

Identification of all possible state variables or those items needed to describe the 
specific tactical situation is the first step. This was accomplished in several ways by ex- 
amining previous and ongoing tactical movement studies, reviewing doctrinal informa- 
tion on movement, and through discussions with groups of selected officers and combat 
modellers familiar with tactical movement. The following 14 variables are identified as 
the significant factors affecting decisions on tactical movement: 

1. Mission 

2. Time available 

3. Equipment and resources available 

4. Threat equipment and capability 

5. Threat intention or mission 

6. Range to the threat 

7. Cover 

8. Concealment 

9. Environment, both weather and obscurrants 

10. Area of operations or theater 

11. Speed or vehicular agility 

12. Range or distance to the objective 

13. Obstacles, both natural and man-made 

14. Artillery 
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The problem now becomes how to organize these variables into scenarios con- 
taining organized combinations or instances in order to provide a means of measuring 
certain attitudes concerning tactical movement. The number of state variables and their 
respective levels must be structured to reduce the required number of responses while 
maintaining the desired degree of fidelity in the experiment. 

One technique is to eliminate any variables not significant to a majority of 
sources. This was the case with the variable artillery since the major impact doctrinally 
concerned obscurrants which fall under environment. In addition, other sources consid- 
ered artillery to influence movement only through the use of improved or special 
munitions which usually fall under obstacles. 

A second technique is to fix selected state variables for a given scenario, espe- 
cially those which are categorical variables and not suited to a regression analysis. The 
respondent or judge is made aware of these in the definition of the scenario and the re- 
sulting scale and function are only valid for that specified combination of fixed variables. 
The following state variables are fixed when using this methodology. 

1. Mission 

2. Threat intention or mission 

3. Area of operations or theater 

4. Equipment and resources available 

5. Threat equipment and capability 

A third technique is to identify which variables may be combined when actually 
considered by the decision maker. This is the case when evaluating cover, concealment, 
environment, and range. All of these combine to determine the degree of inter- visibility 
or "line-of-sight" between the element and the threat. This line-of-sight or LOS and the 
number of lines-of-sight are what become critical to an element manuevering on the 
battlefield. There are also a number of excellent line of sight algorithms available for use 
in combat modelling. 

This same technique is applicable in the case of localized speed or vehicular 
agility. Obstacles, environment, and vegetation combine with the terrain itself to deter- 
mine speed. An element moving over the terrain is not concerned with soil moisture or 
tree stem spacing but with the total effect all these variables have on his speed. The 
combatant evaluates movement based on speed because he knows that the greater his 
speed, the less likely it is he may be engaged and the less time he will require to reach 
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his objective in general. The Condensed Army Mobility Management System. CAMMS, 
provides excellent and highly accurate speed mappings for digitized terrain data. 

The variable “range to threat” also requires further investigation. A combatant 
must evaluate many threats simultaneously and they are often different types with 
varying capabilities. The combatant makes his tactical movement decision based upon 
the greatest known or suspected enemy capability which can influence his mission at that 
moment. This becomes his greatest perceived threat. As mentioned previously, the 
combatant makes a determination of the potential lethality facing him based on the 
types of enemy, their range, and their capabilities based on range and disposition. A 
variation of the range dependency equations used in many Lanchester models will be 
used to calculate the “Relative Effective Range” or RER of each enemy system and al- 
low the analyst to compare the lethality of various systems at different ranges much the 
same as the combatant. These have several mathematical forms studied by Seth Bonder 
and Bob Farrell from 1969 through 1974 [Ref. 13: No page). The minimum RER of all 
threat systems with line-of-sight at a particular location then becomes the range which 
equates to that range given in the survey and used in the model. 

One last observation concerns the variable time. Time actually becomes the de- 
pendent variable for which we wish to determine a relationship. The range to the objec- 
tive is also directly related to and affects the time domain network inherent to GYS. 
These two variables are not included in the set of state variables given for each instance. 

This now reduces the number of state variables to the five fixed variables iden- 
tified previously and the following three key variables. 

1. Localized speed 

2. Number of lines of sight 

3. Range to the greatest threat 

The three key variables must now be broken down into distinct levels in order 
to construct the set of instances needed for the opinion survey. An extract of the actual 
survey used is located at Appendix A. A conscious decision was made to keep the re- 
quired number of responses at a level which allowed the the survey to be completed in 
an hour or less by a respondent. The variables are each divided into four levels which 
cover the majority of conceivable occurrances for each variable. This results in 4 3 or 64 
instances to be evaluated and which must be appropriately organized for the survey. 
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2. Survey Construction and Execution 

This is only a pilot survey given to a limited population of U.S. Army and Ma- 
rine Corps officers assigned to the L'.S. Naval Postgraduate School, TRAC- Monterey, 
and Fort Ord. A total of 62 surveys were distributed during a one week period allowing 
two weeks for a response. A minimum of 25 surveys was needed and it was hoped that 
50 usable surveys could be obtained. A total of 42 surveys were received at the end of 
three weeks of which 40 usable surveys were compiled to provide the data set located 
at Appendix B. 

The actual survey used to gather responses is located at Appendix A. The state 
variables, divided into contiguous levels as listed, are organized into the 64 combinations 
or instances. These are then randomized to alleviate any possible influence on the re- 
spondents or judges due to ordering of the variables as each instance is evaluated. The 
number of levels and resulting number of instances has been kept as small as possible 
while still providing sufficient fidelity in the data in order to demonstrate the validity of 
this method without undue effort in the data collection. The cross-country or localized 
speed is broken down into levels of 0 to 5 KMPH, 6 to 15 KMPH, 16 to 30 KMPH, and 
31 to 45 KMPH. The number or degree of LOS is divided into no lines-of-sight. 1 line- 
of-sight. 2 or 3 lines-of-sight, and 4 or more lines-of-sight. The range is divided into one 
kilometer bands out to 4000 meters. The exact value used for each variable within each 
band was also determined randomly in the case of speed and range using a pseudo- 
normal random number generator with the median value of the variable in that level as 
the mean. This provides further assurance that judges are not influenced by ordering or 
repeating of values. 

The fixed variables are defined in an operations order given as part of the sur- 
vey. The Central Army Group (CEN'TAG) area of Germany is the setting for an Army 
'92 balanced team with Ml Tanks, M2 IFV's, and attachments to conduct a deliberate 
attack against TS0 and T72 Tank, BMP IFV, and BRD.V12 ATV equipped Soviet units 
in defensive positions. Information on aircraft, helicopter, and artillery support is also 
provided. 
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The introduction of the survey also contains detailed explanations of the value 
being rated and the rating categories with examples of instances that would normally fall 
into that category. The five categories are : 

1. No Movement Potential 

2. Marginal Movement Potential 

3. Effective Movement Potential 

4. Good Movement Potential 

5. Excellent Movement Potential 

This rating of the state of a tactical element with respect to possible tactical 
movement based on the combination of variables or the instance is then combined with 
the corresponding ratings provided by all other judges in order to establish a data set 
from which an interval scale can be determined. This comprises the raw frequency matrix 
located at Appendix B, Section A. 

D. RELATING THE COGNITVE AND PHYSICAL CONTINUUMS 

The ability to accurately relate the cognitive and physical aspects of the battlefield 
varies immensely depending on the action being studied, the scenario, and the target 
population. Numerous methods can be applied allowing for complete translation of the 
cognitive continuum into the physical or structuring the resultant scales in order to 
perform a combinatorial optimization. Two specific methods and the underlying theory 
associated with them are examined in this study. 

1. Computing Interval Scales from Categorical Responses 

The theory involved in deriving interval scales using the Law of Categorical 
Judgements is fully explained by Professor Glenn Lindsay (Ref. 12: pp. 6-12] and W.S. 
Torgerson [Ref. 10: pp. 207-227]. However, a review of the major points and underlying 
assumptions is worthwhile. There are 6 key assumptions when using this technique. 

1. Judges' feelings about the scale value of an instance are normally distributed across 
the population with mean n, and variance o ?. 

2. The categories must be a mutually exclusive set of successive intervals which col- 
lectively exhaust the property continuum. 

3. Five categories are required for rating instances. 

4. Judges' feelings about a category upper bound are normally distributed so that for 
category j, the upper bound would have mean n, and variance a]. 

5. This also implies that judges feelings about the distance between an instance value 
and category bound are normally distributed. 

6. Category bounds have the same variances for all j categories. 
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The technique takes the raw frequency matrix for all instances as mentioned and 
first divides these frequencies by the total number of respondents in order to determine 
the relative frequencies for each instance and category. These relative frequencies are 
then summed across the categories for each instance to derive cumulative relative fre- 
quencies which are related to cumulative values under a normal curve. These cumulative 
relative frequencies for category j of instance i become the normal probabilities, p, r used 
in computing the interval scale. The p w are grouped based on which instances have cat- 
egories with similar ratings and all p tJ greater than 0.98 and less than 0.02 are rolled over 
into the adjacent category in order to reduce the effect and possible influence of outlying 
judges. Treating the p tJ as the area under a Normal (0,1) Curve, the Z matrix of normal 
values can be constructed. Row, column, and grand averages are computed from this 
matrix along with an estimate of the standard deviation. The interval scale values of all 
instances are computed as the grand average minus the value of the row average times 
the standard deviation. A linear transform is then used to obtain the correct distance for 
the scale. 

The 10 step method of interval scale development given by Lindsay is shown 

below: 

1. Arrange the raw frequency data in tablular form where the rows are instances and 
the columns represent categories. Columns should be sorted in a rank order from 
least to most favorable as left to right columns respectively. 

2. Compute the cumulative relative frequencies for each row and record these values 
in a new table. This table is referred to as the P array and all values of greater 
than 0.9S or less than 0.02 are combined with the following category eliminating 
extreme cases. This creates an n by (m-k) array where k is the number of columns 
removed. 

3. Treating the /^values as leftward areas under a Normal (0,1) Curve, find the cor- 
responding value of Z from a normal distribution table. Record these values as a 
new table which is the Z = || z tJ || array. 

4. For each instance, i, in the Z array; compute the row average, z, . 

5. For each column or category, j, in the Z array; compute the column average, b r 
Note that b t is the value of category j's upper bound on the interval scale. 

6. Compute the grand average, b, of all Z values. 



7. Compute B = ^ (bj — b) 2 which is the sum of the squared column differences. 



m—k 



m—k 

8. For each row compute A, = £ (z fJ — z,) 2 , the sum of the squared individual differ- 



ences. 



9. For each instance compute 




-j- , an estimate of v a 2 + c . 
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, for all i. 



10. Finally, for each row or instance compute S = b 




The S, values are the interval scale values of the instances which are on the same 
scale as the category bounds, b t . This now enables us to impart distance and a rational 
origin to the scale through anchoring and linear transformations. The above procedure 
is performed on the raw data at Appendix B. The data is segregated into 7 groups cor- 
responding to how instances were rated by like categories. The first and last groups 
cannot be utilized in determing interval scales since they represent extreme instances in 
the survey, have incomplete Z fJ arrays, and only contain one boundry point. All groups 
have incomplete Z v arrays which is why they are segregated. A tactic which enables us 
to continue with the scale development requires performing the scaling procedure on 
each group and then performing a linear transformation between groups with two or 
more common boundries [Ref. 12: pp. 18-20]. Each group of data is then normalized 
and an interval scale computed utilizing the A Programming Lanquage (APL) program 
at Appendix B, Section G. The APL program “NORM’’ was originally developed by 
Paul M. Crawford during his thesis research [Ref. 8: p. 89]. The results of the normal- 
ization and interval scale for each group are located at Appendix B. Section II. 

2. Converting Interval Scales to Usable Forms 

The major requirement at this point is to transform the interval scale values for 
all groups into one contiguous scale by establishing initial boundry values for one group 
and then performing the linear transformations on the remaining groups. This provides 
consistent values for all instances evaluated on one usable scale. Another necessity in 
order to use the scale in models is to convert the interval scale into a ratio scale by es- 
tablishing a rational origin and anchor point or points which intersect on the physical 
and cognitive scales. This allows us to directly map the two continuums together simul- 
taneously. The boundry conditions applied will begin with Group II where a rational 
origin can be applied at the boundry in common with the extreme points shown at Ap- 
pendix B. Section E. A method to anchor the other extreme end of the scale will also 
be discussed. 
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a. The Natural Origin and Anchoring in Computing Interval Scales 

The extreme cases reflecting excellent potential correspond to cases which 
reflect little or no anti-potential. We can reverse the rating categories if desired ( which 
was done in this case ) so that excellent TMP or no anti-potential is the worst rating and 
then perform the normalization procedure and develop the scales. This is done in order 
to correctly place the rational origin. The physical scale of the battlefield based on 
traversal times between cells or nodes reflects only the effects of constraints on the 
physical mobility of a given system. At long ranges when there is no inter- visibility with 
opposing forces and the system can move with great speed, there is no effect due to the 
cognitive considerations of the battlefield. This zero point or rational origin on the cog- 
nitive scale corresponds to the effect any cognitive rating of these extreme cases would 
have and also the respective value of this boundry on the scale. An arbitrary value of 
100 was selected for the upper boundry of Group II. This now allows for the use of the 
APL program "TRANS” shown at Appendix B. Section G. This program was also ori- 
ginally developed by Paul Crawford [REF. 8: p. 90]. 

b. The Linear Transformation 

The linear transformation performed by this program solves a set of simul- 
taneous equations on each group of the form: 

Lower Bound = a + ( /? x Best Column Average) 

Upper Bound = a + ( /? x Worst Column Average) 

where the column averages are already computed and the bounds are input. Once y and 
/? are known, the instance values and other intermediate boundry’ values can be easily 
computed. These other boundry values along with those input are then used to perform 
linear transformations on all other groups with two or more boundries in common. All 
groups can then be combined in order to form one contiguous scale and the value of the 
extreme cases located at the rational origin may be added to the set of scale values of 
the instances. The set of transformed instance values is shown at Appendix B, Section I 
in the last column. The set of instance values along with its corresponding set of vari- 
able values are combined to form the regression matrix shown at Appendix B. Section 
I. 

3. Regression Analysis of the Interval Scale Values 

Multiple regression techniques are used on the data and derived scale values. 
Two different regression packages, Grafstat and SAS. are used applying backward elim- 
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illation procedures. Plots of the independent variables taken two at a time with the in- 
stance values of TMP are done and then surface fitting procedures applied in Grafstat. 
These plots are similar to scatter plots and used in much the same way. Relationships 
between independent variables are of particular interest. The example plots in Appendix 
C. Section A reflect the linear relationships between pairs of variables. The fitted sur- 
faces appear as planes. Due to this result, multiple linear regression is used to include 
product terms of all pairs of independent variables. Appendix C, Section B gives the re- 
sults of backward elimination using Grafstat and Appendix C, Section C gives the cor- 
responding results of backward elimination using SAS. Each section gives the best 4 and 
5 variable models based on the largest R 2 with the best acceptable levels of significance. 

The respective best regression equations for the four variable case and the five 
variable case are: 

TMP = 1 18.24 -0.015(/?.YG)+ U.O(LOS) - 2.0{SPD) + 0.2(LOS x SPD) {1} 

TMP= 112.32 -0.012(/?.Y<7) + 11. 4(LOS) - 1. 1{SPD) 

+ 0.21 (LOS x SPD) - 0.000 15(KA’G x SPD) 

These were the only two equations found acceptable based on the multiple regression 
analysis. 

4. The Psychometric Power Function 

The “Psychometric Power Function” is the selected best regression equation 
from the regression analysis. The scale can be fully determined by taking the set of 
maximum and minimum values for the independent variables and using these to com- 
pute the corresponding maximum and minimum values of TMP from the regression 
equation. The user must be aware that use of this equation or power function is only 
valid for the initial range of values surveyed and used in the regression. The minimum 
value reflects the rational origin or highest TMP which imparts zero additional time to 
the physical scale as discussed previously. The problem now becomes how to anchor the 
other extreme of the scale which reflects the maximum scale value or no TMP. 

5. Translating Values of the Interval Scale to the Physical Scale 

The essential step proposed by this methodology requires translating the cogni- 
tive values obtained using psychometrics into the same scale as the physical traversal 
time scale of the battlefield or to be able to somehow perform a combinatorial opti- 
mization between the two scales. There is a fortunate structure in the data resulting 
from the survey. Referring to Appendix B, Section E; there is only one extreme instance 
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occuring on the end of the scale opposite the rational origin. It is therefore possible to 
prompt the user or survey a selected population to determine the amount of additional 
time in mission completion a combatant would be willing to accept in order to avert such 
a tactical instance. It is only possible to construct such an anchor due to the fact that 
there is only one or are a relatively small number of related instances placed at this ex- 
treme end of the continuum. The scale value of this combination of variables or instance 
is easily calculated. It becomes quite simple to then translate the cognitive scale to the 
time scale by dividing the scale or TMP values by the maximum value computed from 
the power function for the extreme instance which results in a (0,1) scale. Then multiply 
through by the given additional time value from a population mean or as provided by 
an individual respondent. This time scale obtained using the cognitive data can be used 
additively with the physical traversal time scale. Values for each individual user or a 
mean value for a selected population could be preferred depending on the application. 
This translation allows for speed and simplicity when performing the optimization using 
Dijkstra's algorithm which is discussed in Chapter IV. Such an approach is highly pre- 
ferred if at all possible. I Iowever. it is often not possible to relate the cognitive and 
physical scales as in this research. Other optimization techniques which work in combi- 
natorial optimization applications are also discussed in Chapter IV. 



IV. MODEL CONSTRUCTION 



A. GENERAL 

This chapter outlines overall structure of the model and how the methodology al- 
ready discussed is controlled and implemented. There is also some detailed examination 
of algorithms and methodologies; particularly concerning line-of-sight, relative effective 
range, and minimum cost optimization. This study follows a logical order beginning with 
required inputs, continues with algorithms and theory used to process data in the model, 
examines optimization techniques, and concludes with the model output. 

The research version of this model, called the Tactical Optimum Path or TOP 
Model, is located at TRAC; Monterey, California and at the Mobility Systems Division, 
U.S. Geotechnical Laboratory; Vicksburg, Mississippi. The current model contains ap- 
proximately 80.000 lines of code with most routines done in FORTRAN and the user 
interface, graphics, and data management executed in the C Programming Language 
which is used as a shell. More detailed documentation and a draft users manual for the 
model will be published as a joint technical report through the U.S. Army Waterways 
Experiment Station in late 19S9. 

B. TERRAIN DATA BASES 

The digitized terrain data bases used for this model are the identical terrain data 
bases used for the CAM.MS {Version 2.0). This is 100 meter grid explicit terrain data 
translated specifically for this version of the CAM MS [Ref. 7], The model allows for se- 
lection of a geographic region, country, area, and specified “Quad” sheet. This selection 
corresponds to one of the fixed scenario variables, area of operations, identified previ- 
ously. 

A detailed explanation of explicit terrain data and its use can be found in Lecture 
Xotes In High Resolution Combat Modelling by James K. Hartman [Ref. 4j. We will 
therefore dispense with any major discussion of these data. It is important to note that 
many of the terrain parameters provided in the data set are used in other algorithms. 
This is true of the line-of-sight (LOS) algorithms used. Portions of these algorithms have 
been altered to conform to this specific terrain data base. 
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C. REQUIRED INPUTS 

The remaining required inputs primarily relate to the scenario development or fixed 
variables as well. Two of these variables needed as inputs are the appropriate missions 
of each force. Two other variables comprise the equipment mix available to each force. 
These variables are input to the model by the user during initialization of the program. 
The variables listed up to this point encompass all the fixed scenario variables for the 
methodology. It is important to stress that only the psychometric function determined 
for that specific combination of variables or scenario is valid for that combination in the 
model. 

There are presently three other inputs to this model which must be addressed in our 
discussion. These are the element locations for each force, current environment, and the 
parameter related to the extreme scale value selection for translating the cognitive scale 
which was discussed in Chapter III. Section D. The necessity of this last input is evident. 
The locations and environment are necessary in determing path start, intermediate, and 
end points as well as range and line-of-sight results. The environment also affects mo- 
bility calculations and may be played using real time data once the Airland Battlefield 
Environment (ALBE) system is incorporated in the CAM MS. 

The user inputs these data to include season, current weather, light condition, and 
soil-moisture condition for the environment. The user must input initial locations and 
any required subsequent locations depending on the mission as each type of combat el- 
ement is entered. The user is then prompted as the last step prior to initiation of the 
program to provide the time value corresponding to the extreme case or anchor cited 
previously or to select the expected value derived from the sample population. This 
concludes all required inputs necessary to run the model. 

D. THE MOBILITY MODEL, CAMMS (VERSION 2.0) 

The CAMMS has evolved from research, testing, and analysis conducted over the 
past 15 years. The CAMMS was initiated in 1983 in order provide a capability for pre- 
dicting mobility in a field environment in real time [Ref. 7], There are numerous outputs 
available from the CAMMS relating to mobility as listed in Chapter II. This model will 
utilize two specific output mappings and files from the CAMMS Route Analysis Pre- 
diction: the Vehicle Speed Prediction Map and a version of the Vehicle Time Contour 
Map. These mappings provide the localized cross-country speed or vehicular agility 
which is one of the independent variables necessary for the psychometric function and 
the traversal times along each segment between nodes which comprises the set of arc 
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traversal values for the physical continuum. These output values are valid for the speci- 
fied item or weapon system on that terrain under the specified environmental conditions. 
The weather and environment are incorporated in all calculations through the Soil 
Moisture Strength Prediction which is one of the first routines performed by the 
CAM MS [Ref. 7]. A speed value is stored for all vehicle types for each grid cell or node 
to be used in later calculations. The arc traversal times for segments between all adjacent 
grid cells or nodes are also stored for later use. 

The output provided up to this point by the CAM MS are all mobility modelling 
calculations to be used in the model. The remaining areas we will examine pertain to the 
cognitive considerations affecting movement and performing the optimization. 

E. LINE-OF-SIGHT MODELS 

The model also determines if lines-of- sight exist and to what degree or how many in 
order to obtain the next independent variable. LOS. There arc several excellent means 
by which to determine LOS. A detailed explanation of line-of-sight calculations based 
on all types of terrain data can be found in Lecture .Xotes In High Resolution Combat 
Modelling by James K. Hartman [Ref. 4], We will therefore dispense with any major 
discussion of LOS modelling here. 

The LOS algorithms used in this model are derivations of two LOS subroutines, 
IVSCAX and VISTA, developed and used by the BDM Corporation in the Battalion 
Combat Outcome Model (BCOM) and U.S. Army TRAC respectively. The subroutines 
allow for point to point LOS determinations or inter-visibilitv fan determinations de- 
pending on the desired application and which routine is used. The subroutines are 
modified to use the explicit grid terrain data parameters from the CAM MS terrain data 
base. 

The number of LOS for a given grid cell or node is determined by using the inter- 
visibility fan from each known or suspected threat location to obtain a frequency count 
of the total number of lines-of-sight for each cell. FLOS. This value is also one of the 
independent variables stored for use later in the model. 

F. RELATIVE EFFECTIVE RANGE METHODOLOGIES 

The value of range used in the model requires a more detailed discussion than do the 
other variables. One of the assumptions given initially addressed the fact that the forms 
of the range dependency equations studied and proposed by Seth Bonder and Bob 
Farrell are correct [Ref. 13]. Several such forms and their use are expounded by James 
Taylor [Ref. 14 : pp. 93-99]. One similar form was used by Gordon Clark in DYNTACS 
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in determing the “Difficulty", E, from his set of comparative judgements [Ref. 3: p. 98]. 
While there appears to have been a great deal of past interest in such a methodology; 
there has not been any significant or sufficient analysis of how lethality and particularly 
perceived lethality on the part of a combatant is affected by range. Such a study is well 
beyond the scope of this thesis and is not attempted, fhis is the primary reason the 
stated assumption is made. However, an explanation of the equation forms and why they 
are used in this research is necessary. 

f he basic approach concerning range in this model is to assume the combatant has 
sufficient technical knowledge of the systems on the battlefield in order to accurately 
assess their potential effectiveness against other systems. It is possible to use Probability 
of Hit, Probability of Kill) II it, and Total Probability of Kill values from classified 
AMSAA JANUS model scenarios along with the classified effective range data for a 
weapon system type to fit one of the equation forms previously proposed. This becomes 
the perceived as well as actual lethality curve for that system. The systems all have var- 
ious maximum effective ranges with differing curves and equations. The systems are 
equated through the use of a baseline range which is assumed to be the range of the 
battlefield against which the actual range to maximum effective range ratio or difference 
is adjusted. This results in a Relative Effective Range, RER, which reflects the range 
given in the survey and used in the model. 

Some technical data from unclassified reports, periodicals, and the general popu- 
lation's knowledge files were easily obtained. These sources provided some idea as to the 
shape and resulting type of equations the lethality curves for various systems possessed. 
Ranges were obtained from current military periodicals [Ref. 15]. The next step is to find 
a set of fictitious yet properly structured and ordered parameters for the equation forms 
identified from these sources. A number of potential forms could be evaluated. Actual 
ranges and parameters examined in JANUS data files or other sources are not used in 
this research. 
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The following four equations are examined based on discussions with and previous 
analysis performed by the aforementioned combat modellers: 



RER = RG \ -m r" 
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where, 

RER is the Relative Effective Range used in the model 
RG is the baseline range used for this tactical scenario 
RA, is the actual range to system i 
RM, is the maximum effective range of system i 

u. a, andy are the respective shape scale parameters for each respective equation form 

Equation number 3 is a basic Power form of the range dependency equation similar to 
the actual versus maximum range ratio attrition equation used in Lanchester models 
[Ref. 14: p.9S]. Equation number 4 is an Exponential form which works well when fitting 
non-ballistic weapon systems with fairly constant lethality over range such as current 
generation anti-armor missle systems. Equation 5 is a Half-Cosine form which works 
well for systems that tend to have much smaller lethalities in comparison to most other 
systems such as small arms. Equation 6 is a Full-Cosine form applicable for systems 
which may have minimum as well as maximum effective ranges. [Ref. 14: pp. 96-99] 
Sample graphical results of the plots performed is useful at this point. Figure 1 
shows a comparison of the Soviet T72 main Battle Tank versus an AT5 Spandrael 
Anti-Tank Missle. 
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COMPARISON OF LETHALITIES FOR BALLISTIC AND MISSLE WEAPON SYSTEMS 

RANGES ARE IN METERS 
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Figure 1. Comparison of Relative Effective Range for Soviet MBT vs. ATGM 

The tank is represented by a Power equation with a n of 4.0 and RM of 3200 meters and 
the ATGM is represented by an Exponential equation with an a of 0.000S0 and RM of 
5000 meters. The RER is plotted on a logrithmic scale due to the extremely low values 
of RER at shorter ranges. Low RER implies greater lethality. Note the significant dif- 
ference in Relative Effective Range out to approximately 1200 meters. The tank is per- 
ceived to have much greater lethality at these shorter ranges. No major difference is 
evident between 1400 and 2000 meters; but, the lethality of the ATGM remains more 
constant over range and Relative Effective Range changes in favor of the missle from 
2500 meters out to the limit of the battlefield. 

A similar comparison exists for the M1A1 Main Battle Tank and the M2, 3 
Infantry Cavalry Fighting Vehicle where the tank is a ballistic system and the IFV is a 
missle system based on the anti-tank missle system in addition to a 25 millimeter chain 
gun. This comparison is shown in Figure 2. 
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COMPARISON OF LETHALITIES FOR BALLISTIC AND MISSLE WEAPON SYSTEMS 
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Figure 2. Comparison of Relative Effective Range for U.S. MBT vs. IFV 

The tank is represented by a Power equation with a ,u of 4.25 and RM of 3500 meters 
and the IFV is represented by an Exponential equation with an a of 0.00090 and RM 
of 4000 meters. Note the significant difference in Relative Effective Range signifying the 
tank is the greatest perceived threat out to 2000 meters and the IFV equipped with the 
TOW-2 missle is obviously greater past approximately 2800 meters. 

Complete listings of the systems included, the equation forms used, and resulting 
parameters are shown in Tables 1 and 2. Only the Power and Exponential forms are used 
in this analysis and in the model. Table 1 lists the values of the maximum range and 
scale parameter, a , for systems represented using the Exponential form equation and 
Table 2 lists the values of the maximum range and shape parameter, n , for systems re- 
presented using the Power form equation. 
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NOMENCLATURE 


COUNTRY 


PARAMETER, a 


MAXIMUM 

RANGE 


AT4 ATGM 


USSR 


0.00085 


40(»0 m. 


AT5 ATGM 


USSR 


0.00080 


5000 m. 


AT3 ATGM 


USSR 


0.0030 


1500 m. 


BMP-1 IFV 


USSR 


0.00275 


2000 m. 


BMP-2 IFV 


USSR 


0.00080 


4000 m. 


BRDM-2 ATV 


USSR 


0.00075 


5000 m. 


HIP HELO 


USSR 


0.00180 


4000 m. 


HIND HELO 


USSR 


0.0010 


4000 m. 


MOO 1 1TV 


USA 


0.0010 


4000 m. 


M2 3 IFV CFV 


USA 


0.00090 


4000 m. 


TOW ATGM 


USA 


0.00125 


4000 m. 


DRAGON ATGM 


USA 


0.00250 


1000 m. 


AH IS HELO 


USA 


0.0010 


4000 m. 


APACHE HE I O 


USA 


0.00075 


4000 in. 



Table 1. TECHNICAL DATA FOR MISSLE SYSTEM EQUATIONS 



It is also important to note that we can compare NATO and Soviet systems using 
these curves. Note that the T72 Main Battle Tank has a lower ,u value and is not as 
effective as the M1A1 Main Battle Tank nor does it possesss as long a maximum effec- 
tive range. The difference is noticeable when the plots shown in Figures 1 and 2 are 
compared. 

It is also interesting to note how the lethality for various systems change with re- 
spect to range. This gives some idea as to whether or not the RER calculated makes 
sense or is verifiable. The small change in slope of approximately 0.1 observed within 
the first 1000 meters of actual range for both tank weapons systems is one example. This 
range is within the boresight range of each system and there is little or no change in the 
trajectory of tank cannon munitions at such short ranges. The slope and subsequent 
RER support these facts. The slope between 2000 and 3000 meters actual range in- 
creases to approximately 1.7 indicating perceived lethality decreases more quickly. There 
is almost no difference in the RER of 1 at 500 meters range and the RER of 14 at 1000 
meters range for the M1A1 Tank. However, there is a significant difference in the RER 
of 2600 at 3500 meters for the tank. The TOW-2 missle on the IFV has an RER of only 
1200 at 3500 meters range. While these equations and parameters are fictitious, their 
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application and results correspond to the manner in which experienced military person- 
nel tend to compare these systems. 



NOMENCLATURE 


COUNTRY 


PARAMETER, n 


MAXIMUM 

RANGE 


T54 55 MBT 


USSR 


3.50 


2000 m . 


T62 64 A MBT 


USSR 


3.75 


2S00 m. 


T64B 72 MBT 


USSR 


4.00 


3200 m. 


TSO MBT 


USSR 


4.25 


3200 m. 


BTR60 60 P PC 


USSR 


0.50 


1500 m. 


DIS. INF. 


USSR 


0.40 


1000 m. 


UNK. VEH. 


USSR 


1.50 


4000 m. 


UNK. OTHER 


USSR 


1.00 


4000 m. 


M48 60 MBT 


USA 


3.80 


2S00 m. 


M 60 A 1 MBT 


USA 


3.90 


3200 m. 


M60A3 MBT 


USA 


4.00 


3200 m. 


M1P1 MBT 


USA 


4.10 


3200 m. 


MIA1 MBT 


USA 


4.25 


3500 m. 


Ml 13 PC 


USA 


0.50 


1 200 m. 


M60CEV 


USA 


o.so 


1000 m. 


DIS. INF. 


USA 


0.40 


1000 m. 


UNK. VEIL 


USA 


1.50 


4000 m. 


UNK. OTHER 


USA 


1.00 


4000 m. 



Table 2. TECHNICAL DATA FOR BALLISTIC SYSTEM EQUATIONS 



For this model, an RER value for each cell or node is determined by taking the ac- 
tual range, RA, for each system with line-of-sight to that cell and using this range in the 
appropriate RER equation form for that system. The minimum value of the RER for 
all systems with line-of-sight, if more than one exists, is retained as the RER value for 
a cell. 

The data and equation forms identified and derived from this short study do need 
to be evaluated more thoroughly. However, such an effort is beyond the scope of this 
thesis and would require development of an extensive experiment in order to sufficiently 
determine how perceived lethality is affected by range for various weapons systems. 
Equation forms with exact parameters derived from existing AMSAA data would almost 
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certainly be classified as would the results of such an experiment. This is certainly a 
proposal for future study and enhancements: but. the assumption stated is used at this 
point in order to continue development of methodology based on the use of 
psychometric techniques without undue requirements due to an analysis of perceived 
lethality and or inclusion of classified data. 

All required inputs to the model are complete once the CAM MS has been run; the 
LOS calculations are performed; the SPD, FLOS, and RER values have been obtained 
for each cell; and the arc traversal times between all adjacent cells are stored. The ap- 
plication of the cognitive portion of the model and transformation of the values to be 
used in the optimization can now be done. 

G. APPLYING THE PSYCHOMETRIC FUNCTION 

The physical time scale is already determined from the CAM MS but we must utilize 
the other data as inputs to the psychometric or decision function derived from the re- 
gression analysis in order to obtain the cognitive or tactical time scale. The regression 
equation or decision function uses FLOS, RER, and SPD as independent variables in 
order to output a dependent variable, Tactical Movement Potential or TMP. This value 
is computed as a ratio scale value for each cell through application of a rational origin 
and then translated into the corresponding physical time scale made up of the traversal 
times through the use of an anchor as outlined in Chapter 111. Section I). We now have 
two sets of corresponding time values. 

The cognitive or tactical time values represent the additional cost in units of time, 
equivalent to those time units on the physical scale, for a combat element to enter that 
cell based on the existing set or combination of state variables. The two sets of time 
values are additive. Consequently, the tactical time computed for each cell is then added 
to the physical traversal times stored for each arc which enters or has a tail node in that 
particular cell. This set of time values becomes the total time along each arc and the set 
of values desired for use in the optimization. 

H. THE OPTIMIZATION 

Optimization of this set of translated network values is performed using Dijkstra's al- 
gorithm which is relatively quick and efficient in this application. However, it is not al- 
ways possible to determine a rational origin and particularly one or more anchors in 
order to translate the cognitive scale. Dijkstra's algorithm can solve the optimization 
problem very simply if the requirement is to either reach the objective by minimizing 
tactical difficulty which is equivalent to maximizing TMP. to reach the objective by 
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minimizing only the traversal time, or to reach the objective by minimizing some total 
translated time as in this case. However, the problem is one of combinatorial opti- 
mization when it is not possible to determine a relationship between the two scales. We 
will first conduct a quick review of minimum cost optimization using Dijkstra's algo- 
rithm and then examine one method for performing the same optimization combina- 
torially on separate physical and cognitive scales. 

1. Dijkstra's Algorithm 

Dijkstra's algorithm uses a label setting procedure on a network graph with 
non-negative costs for all arcs in the network. This network structure must be main- 
tained in order to avoid possible negative cycle lengths. Two sets of cells are maintained. 
The first set are those cells whose shortest distance or minimum cost path is already 
known. The second set consists of all other non-traversed cells. The known set is initially 
comprised of only the start or source node. At each iteration, a new cell is added to the 
known set which has the shortest distance or minimum cost path to the start point. The 
minimum cost value or shortest distance can be recorded or a predecessor array main- 
tained. Eventually, all possible cells on the battlefield are included in the known set. At 
this point, the minimum cost path from the source to any other cell is available and the 
value can be determined. [Ref. 16: pp. 203-20S] 

Dijkstra's algorithm can be made even more efficient by limiting the nodes 
searched at each iteration. This is applicable in the model since the only permissible arcs 
are those between adjacent cells. The algorithm can also be used to perform a minimum 
cost optimization along a set of intermediate objective points or cells through applica- 
tion of the “Sliding Pattern Concept” [Ref. 5: pp. 22-26], The optimization is performed 
locally between the ordered set of cells n times, where n is the number of intermediate 
cells when the order is specified and must be maintained. The optimization is performed 
in total until all intermediate cells are labeled or in the known set and all are contained 
along one unique path. This can be accomplished using any number of data structures. 
The set of objective points or cells is input by the user when the model is being initial- 
ized. In summation, Dijkstra's algorithm is quite useful for almost any network and can 
be adapted to make it even faster and more efficient with little work. 

2. Combinatorial Optimization 

The optimization problem can be solved very simply if the requirement is to ei- 
ther reach the objective by minimizing some cost along one unique scale when the cog- 
nitive and physical scales can be related as demonstrated in the previous sections. 
However, the problem is one of combinatorial optimization when it is not possible to 
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determine a relationship between the two scales. This requires the analyst to find a 
method by which it is possible to minimize tactical difficulty while also optimizing the 
traversal time. The optimized time must be less than time T and the tactical difficulty 
must also be minimized within this optimization. Time T should be given by the user 
or obtained from the mission statement of the scenario. The model should select a route 
which satifies all constraints while solving one equation for both tactical difficulty, td x , 
and the traversal time, t x . A method that combines both constraints into one equation 
is the use of a Lagrangian multiplier, A > 0, which can be used in the following equation: 

min(f</ +/x t)x— A x T {7} 

all x 



where. 

x is a vector of the arcs which make up a network solution 

td is the tactical difficulty based on exposure associated with each .v 

t is the travel time associated with each .r 

/ is the Lagrangrain multiplier 

T is the maximum time limit to traverse the network 

a. Assumptions 

1. T is larger than the time it takes an element to reach its objective using a shortest 
time network solution. 

2. The network is undirected. This allows an element at destination B to backtrack 
along its route if that route offers a less difficult route and still meets the T time 
limit to reach the final objective. 

b. Methodology 

The equation is used to calculate the shortest path based on A using 
Dijkstra's algorithm. The algorithm shown at Appendix I uses a hierarchial adjacency 
list to search for the shortest path between the start point and the intermediate point 
and then from the intermediate point to the final objective. Other data structures could 
be used. 

To calculate the maximum possible value, ). mw , the largest t x from all the 
edges in the network is multiplied by the number of vertices in the network. If A m „ is 
used as ). in the equation, the objective function minimizes the time traveled throughout 
the network with no regard to tactical difficulty. When A is zero, the minimum possible 
value of tactical difficulty is found. There exists a A such that 0 < A < A miX gives the op- 
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timum tradeoff in minimizing the tactical difficulty while ensuring the traversal time 
constraint is not violated. 

The best lower bound for the equation is determined by solving for 
maxL(/) which can be solved bv doing a binarv search from / min to / max . As long as the 
slope of L(X) is non-positive then the route selected is a feasible solution with regard to 
T. A binary search is conducted on the interval between / mm and / max until a specified 
interval, e . is obtained between / max and Z m , n ., This X with e sufficiently small indicates 
the search for the optimum path minimizing tactical difficulty while arriving at the ob- 
jective within T traversal time has been accomplished. 

c. The Combinatorial Optimization Algorithm 

A VS-Fortran program was developed to solve the tactical movement 
problem (see Appendix I) and consists of the following steps: 

1. The data set containing the edges and their associated time and difficulty values are 
read into a hierarchial adjacency list. The value of the maximum t x and TDM AX 
in the network are found. 

2. Inititially Dijkstra's algorithm is called to find the slope, objective function value 
and path for / = 0 and / max = TDM AX x (Xocles— 1) 

3. A check is conducted between z mm and / max and if X m3K — X mm < e then the program 
is completed and the shortest path with the lowest tactical difficulty has been 
found. 

4. X,„, d is the midpoint between z min and z max 

5. Dijkstra's algorithm is called to compute a slope, objective value and path based 
on X md 

6. If the slope is positive then X mid becomes z min and Step 3 is repeated. Otherwise the 
slope is non-positive, indicating a feasible path and X mi<t becomes the new z max . The 
objective function value of X md is compared to the previous best objective function 
value, OBJMAX. If ). md objective value is larger than OBJMAX, ). md objective 
value becomes the new OBJMAX and the optimal path, OPTX, is updated. Go to 
Step 3. 



Once the e value has been reached the optimum path and associated op- 
timal values for time and tactical difficulty are outputted to a file. 

</. Conclusions 

The complexity of Dijkstra's algorithm is 0( | V\ 2 ) and the complexity of 
the binary search is 0{\og 2 TDMAX). The Dijkstra algorithm is utilized with each iter- 
ation of the binary’ search; therefore, the total complexity of the program is 
0( | V | 2 \og,TDMAX) 
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equation 



To calculate the required number of iterations for a set of nodes, use the 
to determine the interval required and then solve for the number of 
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iterations, k, using the equation 



k > 
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or use an e of sufficiently small value. 

Changes to the program can be made to decrease the run time if more in- 
formation is known about the nature of the network. In a directed network the 
Dijkstra's algorithm can be changed to stop the search once point B has been found and 
then continue from point B until point C is found. The search could also be restricted 
over certain portions of the network based on the network's structure, greatly reducing 
the number of searches. Another method that could be incorporated is changing the 
binary search to a slope intercept search over the interval between / min and ). mn . These 
methods were not incorporated in this model since it could be used to model several 
types of networks with varying numbers of intermediate points and structures. 

I. MODEL OUTPUT 

The output for the model is designed to be presented graphically in addition to the 
data files for use in models. Some output is available interactively during model initial- 
ization. All output is available upon completion of a model run for use in planning, 
analysis, and to provide an audit trail. The key output is certainly the set of optimum 
path points and total route traversal time. 

Interactive output consists of LOS fans displayed to the screen, a selection of base 
maps which include terrain elevation or vehicle speed maps, element locations, area of 
operation outline, and path point locations. Upon initiation of the model, cells are la- 
belled in accordance with Dijkstra's algorithm and shown on the screen until completion 
of the routine, at which time the optimum path is displayed. The data are also stored for 
later use. Elimination of the interactive graphics capability when the model is interfaced 
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with a combat model will greatly reduce run time. 

Final output is comprised of a selection of mappings which are the vehicle speed 
map or SPD map, FLOS map. RFR map, TMP map, and total time map. The locations, 
path points, area of operations outline, and actual optimum path can be imposed on any 
one of the selected base maps. Base maps for hard copy output may be selected in scales 
from 1 to 25,000 up to and including 1 to 250,000. Actual cell values for each state var- 
iable can be accessed by selecting a point interactively with a bit mapping utility. The 
set of optimum path points are easily obtained from a data file. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. GENERAL 

The working model is resident in the warlab of The TRADOC Analysis Command; 
Monterey. California. Initial verification runs were conducted on 17 through 25 July 
1989. The model is specifically tailored so as to optimize speed and efficiency on a 
PC-AT 3S6 computer. A Dell System 325 with VGA-Plus graphics and a high speed, 
wide carriage color printer for hard copy output was purchased to support the model. 
Results and subsequent conclusions and recommendations based on the methodology 
implementation and model runs are addressed in the following two sections. 

B. CONCLUSIONS 

I. Model Verification 

The ability to perform a validation of the model within the scope of this thesis 
is impractical if validation is possible at all. However, verification in terms of the general 
behavior exhibited by the model and the sensibility of the results must certainly be ex- 
amined. A very strong verification is not even possible without additional support and 
experimentation. Certainly, a more stringent verification or validation should be the 
focus of additional study; but. we must rely on the knowledge of experienced military 
officers and combat modellers in order to accomplish this limited verification. 

There were 16 different variations of fixed state variables initiated for the test 
scenario, each with a different set of beginning path points and objectives. The expected 
value of 5.5 minutes for the scale time value is used in 10 of the model runs. A value of 
0.0 minutes is used in 2 of the model runs in order to compare the results to the base 
CAMMS model runs where there is no effect due to TMP. There are no discrepancies 
in these comparisons. There are 4 other model runs which use times varying from 2.0 to 
100.0 minutes in order to determine some idea of the sensitivity of the model with respect 
to the scale time value. This is related to allowing the user to input this value based on 
their own risk posture where a large value implies risk aversity. Some output from these 
verification runs are located at Appendix E. A full scale demonstration of the model was 
presented to a selected audience of military officers and professors on 19 July 19S9. 
Impressions concerning the model results and and proposed enhancements were solicited 
with an enthusiastic response to the results. 
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2. Output Analysis 

The model appears to work extremely well in terms of the route selection proc- 
ess using the expected value and in terms of the sensitivity analysis for varying risk 
postures. The first four pairs of 1 to 25.000 scale plots shown in Appendix E reflect 
movement based strictly on mobility versus tactical movement using the expected scale 
time value. The last two pairs reflect a change in risk posture where the expected time 
scale value is compared to greater times signifying a more risk averse posture. 

The results are quite evident from these output mappings. The methodology and 
use of a psychometric decision function does provide a set of path points which appear 
to optimize TMP during movement in addition to mobility considerations. The model 
can also be applied deterministically, stochastically, or individually using different values 
for the scale time value. These output also appear to support sensitivity of the model to 
varying risk postures. 

C. RECOMMENDATIONS 

1. Enhancements and Areas for Research 

a. Approved Scenario List 

The scenarios discussed in the development of this methodology are a key 
element in further applications and research concerning the methodology. A complete 
set of scenarios of interest detailing all possible combinations of the fixed state variables 
must be determined in order to structure a data collection effort and target appropriate 
sample populations. There will eventually be a decision function derived for each of the 
identified combinations if the methodology is to be used. This effort would be best per- 
formed through the elforts of the U.S. Army Center for Army Lessons Learned (CALL) 
using information from throughout the military. 

b. Expanded Survey and Regression Capability 

The fidelity of this survey and corresponding regression were discussed pre- 
viously. Continuing data collection and analysis elforts should be more detailed in terms 
of categorization and the number of levels or groups used for the independent variables. 
There is certainly a balance needed between the amount of fidelity desired in the model 
and the magnitude of the collection and analysis effort. The populations used for these 
efforts should also be expanded and examined closely to insure the appropriate sample 
is surveyed. Larger samples would certainly lower the error and increase the power of 
any such test. These actions should be a cooperative effort between the Army Research 
Institute (ARI) and TRADOC Analysis Command (TRAC). 
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c. Multiple Vehicle Paths and Formation Structure 

The need to model units using a formation structure or by defining a path 
width for each arc on the battlefield which can be used to control movement will be a 
necessary requirement prior to interfacing this model with other combat models. This 
area of research alone is practically unbounded. A decision should be made by the Army 
Model Board concerning which approach to use. A formation structure could be applied 
similar to that found in DYN'TACS [Ref. 3; pp. 150-155] or developing path parameters 
for the digitized terrain data and then limiting the optimization to only those arcs with 
a minimum acceptable width [Ref. 17: pp. 33-36] are both feasible. It is also possible to 
perform a combinatorial optimization using LaGrangian relaxation on three variables 
to include path width. This effort would be best accomplished by the T RADOC Analysis 
Command (TRAC). 

</. Relative Effective Range Determinations 

A similar type of pyschophysical analysis needs to be conducted in order to 
determine if combatants' perceived lethality sufficiently corresponds to actual lethality 
curves. The same method of categorical judgements could be used to derive exact re- 
gression equations or lethality curves for desired weapons systems in the model or actual 
classified lethality curves could be instituted if there is strong correlation between the 
perceived and actual results. The data and curves would certainly be classified in either 
case. These actions should also be a cooperative effort between the ARI and TRAC or 
the Army Material Systems Analysis Activity (AMSAA) and TRAC. 
e. Symmetric Applications 

There is also a need to determine if psychometric analysis and functions 
done with U.S. and allied sample populations can be extrapolated to and or directly 
applied to Soviet or other forces and if so. at what levels and for which actions. This 
analysis is necessary unless some drastic assumptions are made concerning symmetry 
between opposing force decision processes and tactics in general. This effort will cer- 
tainly be classified and require the assistance of all major intelligence agencies. 

2. Methodology Utilization and Future Applications 
a. Tactical Decision Tool 

A major advantage of this methodology used in a model versus a Tactical 
Decision Aid (TDA) such as the CAMMS is the ability of this model to provide or make 
a tactical decision in real time. The CAMMS and other models only provide pertinent 
information to the combat decision maker to assist him in choosing a course of action. 
The model can make tactical decisions in a laboratory or analytical environment and 



40 



could be used in the field environment given the appropriate terrain data bases and in- 
telligence input. The system and model have been developed to be easily transportable. 

b. Tactical Training Tool 

The model would make an excellent Tactical Training Tool allowing per- 
sonnel to compare their decisions in a laboratory or classroom environment to computer 
simulations based on the model. This application could be performed using a variety of 
mediums. This application could even be done as a field exercise given the flexibility of 
current personal computers. 

c. Combat Model Preprocessor / Subroutine 

The methodology could be used as a pre-processor to provide sets of opti- 
mum path points during the initialization of any number of current combat models and 
then be used interactively as an event driven subroutine to recompute path changes 
during the course of a battle. This is highly desirable in development of systemic models 
where we wish to remove the human from the process and decrease variability in the 
battle calculus. The major problems will be to determine which events would drive an 
update in the path evaluation and whether to use the paths at an individual vehicle or 
unit level depending on resolution. 

The model has been designed and structured to make it compatible with 
several existing combat models in terms of the terrain data bases used, the path output, 
and program lanquages used. These include the Battalion Combat Outcome Model 
(BCOM) developed by the BDM Corporation, both versions of JANUS combat models. 
CASTFOREM, and others. 

d. Research and Analysis Tool 

The methodology has numerous applications and great potential as an 
analysis tool in current form or as part of a systemic research model. The methodology 
can be applied in evaluating any type of decision process when the appropriate state 
variables can be determined and a proper sample population identified. Similar tech- 
niques can be used to include these decision functions in models leading to true systemic 
models and expert systems. 

This specific model, once scenarios are determined and analysis performed 
to fully expand the model, can be used to analyze alternative General Defense Plans 
concerning movement. It may be incorporated into models such as the Obstacle Planner 
System (OPS) to assist in analyzing optimum placement of available engineer resources 
in countermobility or mobility roles. These are only a couple of possible research and 
analysis applications. 
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e. The Generalized Value System, GVS 

The methodology also provides a theoretically sound method of determining 
total traversal times for the time domain network of the battlefield which is a necessary 
input to the Generalized Value System. There is additional work to be done in applying 
the methodology to the GVS. This specifically relates to the decision level at which a 
future state decision is made in the GVS versus the level at which paths are computed 
in the model. The problem does not exist at lower levels such as section or platoon. 
However, it may not be feasible to extrapolate values of the time domain network at a 
higher level from a number of smaller units comprising a large unit. The use of the path 
width discussed previously may be one solution. However, there would need to be a 
method of determining traversal speeds and other data for aggregated units with com- 
binations of vehicles which the CA.MMS and this model can not accomplish. 

D. SUMMARY 

The methodology demonstrates a way to determine optimum, tactical movement 
paths for a specified vehicle and. or small unit based on the operator's cognitive decision 
processes, as well as the physical effects of terrain and environment on mobility. The 
approach uses psychometric techniques inherent to the Generalized Value System (GVS) 
in order to determine a decision function based on the specific tactical scenario and given 
equipment configuration which provides a means to determine the Tactical Movement 
Potential (TMP) for each terrain cell. This cognitive value in an interval scale can then 
be translated into the same scale as the physical continuum. The cognitive time value 
based on the user's decision process is then added to the physical traversal times for each 
cell computed from output provided by The Condensed Army Mobility Management 
System (CA.MMS). This renders a value mapping which can be optimized using 
Dijkstra's Algorithm or a combinatorial method could be instituted. The resulting sets 
of path points are optimized for speed time and the cognitive tactical considerations 
evaluated using these psychometric methods thereby rendering a truly optimum, tactical 
movement path and optimum traversal time. The movement path and resulting times 
can be used in a multitude of modelling or analytical applications and particularly in 
determining the time values needed to compute the Situationally Inherent Power (SIP) 
of the GVS. In actuality, this methodology could be applied to almost any tactical de- 
cision process in the development of expert systems and models and shows great poten- 
tial for future use. 
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APPENDIX A. A SURVEY OF FACTORS INFLUENCING TACTICAL / 

COMBAT MOVEMENT. 



A. INTRODUCTION. 

1. Purpose. 

The purpose of this survey is to obtain an estimate of the degradation to a sys- 
tem's Tactical Movement Potential (TMP) based on changes in 3 key variables and 
combinations of these variables. The results of analysis to be performed on these survey 
responses will then be used to develop input to a larger tactical decision tool designed 
to provide sets of optimum path points to commanders in the field. Your response to 
each question will reflect the value you plaee on a given combination of the 3 key vari- 
ables with respect to Tactical Movement Potential (TMP). Your responses as a tactical 
decision maker commander will assist in determining accurate representations of how 
changes in each key variable affect the rating of its relative importance to tactical 
movement. They may also enable us to determine how the key variables are correlated. 

2. Structure. 

a. Key Variables. 

There are 3 key variables which have been identified as critical in determin- 
ing a system's TMP. These 3 key variables are: 

1. System Speed based on the cross-country movement rate. 

2. Cover and Concealment based on enemy line-of-sight (LOS). 

3. Range based on the distance to known or suspected enemy positions. 

b. Rating Categories. 

There are 5 rating categories from which to choose in evaluating TMP. 
These are a mutually exclusive set of successive intervals that collectively exhaust all 
possible values of the TMP. These 5 rating categories are: 

1. No Tactical Movement Potential. 

2. Marginal Tactical Movement Potential. 

3. Effective Tactical Movement Potential. 

4. Good Tactical Movement Potential. 

5. Excellent Tactical Movement Potential. 
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c. Discussion. 



A short discussion to insure complete understanding of the variables and 
the rating category you should select based on your tactical experience and judgement 
regarding each scenario is required. We will first examine the key variables and then 
define discuss the rating categories. 

The 3 key variables represent the primary factors considered by a vehicle 
or unit commander when selecting an Axis-of-Advance or Route-of-.Movement in oper- 
ations such as an Attack, Movement-to-Contact, Withdrawal, or other operation. In 
each given scenario; which is simply a combination of the 3 key variables; the system 
employed will have a varying TMP. For example, a system travelling over rolling, open 
terrain at a range exceeding 4 km. from the nearest opposing force with a formation of 
low hills preventing enemy observation based on LOS would be capable of moving at 
great speed and in relative safety. This would result in high values for each key variable 
and exhibit "Excellent Tactical Movement Potential". On the other hand, a system 
moving through a peat bog and becoming virtually mired within 800 m. of several known 
enemy positions with little or no cover and concealment would almost certainly be de- 
stroyed or at least badly damaged and become combat ineffective. This would result in 
low values for each key variable and suggest the system has "No Tactical Movement 
Potential". Varying combinations of the 3 key variables render scenarios which will be 
rated between these two extremes. This leads to the second portion of our discussion. 

In order to insure consistent responses between the respondents when 
evaluating each scenario, all respondents must thoroughly understand the definition of 
each rating category. These categories arc: 

• "No Tactical Movement Potential" implies the system cannot move through an 
area, a No-Go area such as a swamp, or that the system has practically no proba- 
bility of surviving such a scenario in the judgement of the respondent. This may 
be due to a slow speed, numerous enemy observations at a close range, etc.. 1 his 
represents a tactical situation any prudent commander would never select if possi- 
ble. 

• "Marginal Tactical Movement Potential" implies the system can move at a mar- 
ginal rate and has some probability of surviving the scenario which is low. How- 
ever, this is a situation where although the system can survive; the simple act of 
survival would render the system combat ineffective for a period of time. This re- 
presents a tactical situation which is only acceptable as a last resort for brief peri- 
ods of time. 

• "Effective Tactical Movement Potential" implies the system can move reasonably 
well and or has a moderate probability of survival and mission success all things 
considered. This represents a tactical situation most commanders would consider 
acceptable in most cases although not preferable. 
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• "Good Tactical Movement Potential" implies the system can move at better than 
average cross-country movement rates and or has a relatively high probability of 
survival and subsequent mission success. This is certainly a preferable tactical sit- 
uation. 

• "Excellent Tactical Movement Potential" implies the system can move at or near 
the maximum cross-country speed with relatively little or no probability of being 
engaged. This is the optimum tactical situation. 

3. Instructions. 

The survey consists of three sections. This 'Introduction' is the first section fol- 
lowed by a 'Personal History ' section used to collect population data. The actual tac- 
tical situation and survey compose the 'Tactical Movement Survey' section which is the 
last but most important section. 

Please insure you feel comfortable with the variables and especially the rating 
categories. Fill in the required personal data as completely as possible and proceed to 
the survey. 

You will begin by reading a synopsized OPORD giving the Task Organization, 
Situation, and Mission. Once you fully understand the tactical situation you will con- 
tinue by answering or rating the TMP for 64 questions scenarios. Please rate ALL 
scenarios in accordance with your initial judgement in as timely a manner as possible. 
This is done in order to replicate the manner in which such decisions are made in fast- 
paced. tactical operations. Do not change answers once you have recorded your initial 
response unless the change is due to an inadvertent error. Answer ALL scenarios. Thank 
you for your cooperation. 
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B. PERSONAL HISTORY DATA. 



1. Present Rank: 

• LTC, 0-5 or above 

• MAJ, 0-4 

• CPT, 0-3 

• LT, 0-1 & 0-2 

• E8 - E9 

• E6 - E7 

• E5 or below 

2. Branch or MOS: 

• Armor Cavalry 

• Infantry 

• Aviation 

• Air Defense 

• Field Artillery 

• Engineer 

• Other 

3. Experience: 

• Time on active duty Years Months 

• Time spent in Combat Arms units Years Months 

• Time attached to Combat Arms units ... Years Months 

4. Opinion: 

Do you feel any other variables should be included in the scenarios which could 
significantly affect the tactical decision logic? 

• Yes No 

If so, what are these variables and how do they impact tactical movement? 
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C. TACTICAL MOVEMENT SURVEY. 

1. OPORD 1-89, HQ/TF 1-15 XXX 

Task Organization : 

Team Yankee 1-15 XXX (-) 

2 Mech. Pit. (M2) 2 Tank Pit. (Ml) Team HQ's 

1 Engineer Sqd. w, CEV 1 Stinger 1 Fist 

1. Situation. 

a. Enemy Forces: 

Unidentified forces of the Soviet 111 TIT MRR are defending in sector. 
The enemy is estimated to be at 72 % strength in men. equipment, and supplies. 
They are preparing defensive positions. The enemy is equipped with organic 
BMP IFV's and BRDM ATY's as well as being reinforced by a TSO Guards 
Tank Battalion and what is left of the Regimental T72 Tank Battalion. Expect 
enemy artillery support from the RAG using 122 SP's. Other artillery support 
is possible. The enemy may receive Helo support in the form of IMP or 
HIND-D units if our operations are successful. The status of enemy fixed-wing 
support is not known at this time. The enemy has positions about 5 Kilometers 
away. Some enemy forces have been located and identified (assume you have 
an Intelligence Overlay). Other larger enemy forces are known to be preparing 
positions behind their security zone. 

b. Friendly Forces: 

Team Yankee has been conducting a Movement-to-Contact as part of 
TE 1-15 XXX and the 10TH Brigade. The Scouts and Team X-Ray are to our 
front and have encountered enemy forces in a security zone forward of their 
main defensive belt. Team X-Ray lost 3 vehicles today. Team Whiskey and 
’I earn Zulu are on arc left and right respectively. The TE has a Battalion of 155 
SP's in DS. Teams Whiskey and Zulu will attack abreast to secure intermediate 
Objectives Erick and Frack. 'leant X-Ray will provide supporting fires from 
their present positions. 

2. Mission. 

Team Yankee TE 1-15 XXX will attack XI T 1000 hours passing through 
Team X-Ray vicinity CP66 and between teams Whiskey and Zulu vicinity CP99 in 
order to sieze high ground vicinity Objective Sam. Team Yankee will then defend 
in place until relieved. 



4 . * 



Based on the above given tactical situation; place yourself in the position of a 
vehicle, platoon, or unit commander and give your rating of the Tactical Movement 
Potential for each of the following scenarios. 
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2. Example Survey Scenarios for Questions 34 through 41. 



• 34. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 

» Range to the nearest known or suspected enemy position is about 1475 meters. 

■ Your cross-country speed is less than 5 KMP11 or 3 MPH. 

■ You are within the LOS of 1 known or suspected enemy location. 

What do you rate the Tactical Movement Potential of this scenario as - 

(JNo TMP (JMarginal TMP (JEfTective TMP (JGood TMP (JExccllent TMP 

• 35. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 

• Range to the nearest known or suspected enemy position is 3740 meters 
» Your cross-country speed is greater than 45 KMPH or 27 MPH. 

■ You are within the LOS of 4 or more known or suspected enemy locations. 

What do you rate the Tactical Movement Potential of this scenario as - 

<_)No TMP (JMarginal TMP (JEfTective TMP (JGood TMP ( JExccllent TMP 

• 36. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 

« Range to the nearest known or suspected enemy position is about 750 to 770 
meters. 

• Your cross-country speed is approximately 20 KMPH or 12.5 MPH. 

» You are not within the LOS of any known or suspected enemy locations. 

What do you rate the Tactical Movement Potential of this scenario as - 

(JNo TMP (_) Marginal TMP (JEfTective TMP (JGood TMP (JExcellent TMP 

• 37. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 

» Range to the nearest known or suspected enemy position is 2575 meters 
» Your cross-country speed is approximately 42 KMPH or 26 MPH. 

• 't'ou are not within the LOS of any known or suspected enemy locations. 

What do you rate the Tactical Movement Potential of this scenario as - 

(JNo TMP (JMarginal TMP (JEfTective TMP (JGood TMP (JExcellent TMP 
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• 38. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 

• Range to the nearest known or suspected enemy position is 1440 to 1515 meters. 

■ Your cross-country speed is less than 5 KYI PH or 3 Y1PI1. 

■ You are within the LOS of 4 or more known or suspected enemy locations. 

What do you rate the Tactical Movement Potential of this scenario as - 

(_)No TMP (JMarginal TMP (JEfiective TMP (JGood TY1P (JExcellent TMP 

• 39. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 

• Range to the nearest known or suspected enemy position just more than 3500 
meters. 

• Your cross-country speed is less than 5 KYIPI1 or 3 MPH. 

■ You are within the LOS of 1 known or suspected enemy location. 

What do you rate the Tactical Movement Potential of this scenario as - 

(JNo TMP (JMarginal TMP (JEfiective TMP (JGood TMP (JExcellent IMP 

• 40. You arc moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 

• Range to the nearest known or suspected enemy position is slightly less than 
650 meters. 

■ Your cross-country speed is less than 5 KMPL1 or 3 MPH. 

■ You are within the LOS of 4 or more known or suspected enemy locations. 

What do you rate the Tactical Movement Potential of this scenario as - 

(JNo TMP (JMarginal TMP (JEfiective TMP (JGood TMP (JExcellent TMP 

• 41. You are moving cross-country as a part of Team Yankee attacking to seize 
Objective Sam under the following conditions: 

» Range to the nearest known or suspected enemy position is between 22S0 and 
2320 meters. 

■ Your cross-country speed is approximately 40 KMPH or 25 MPH. 

• You are within the LOS of 4 or more known or suspected enemy locations. 

What do you rate the Tactical Movement Potential of this scenario as - 

(JNo TMP (JMarginal TMP (JEfiective TMP (JGood TMP (JExcellent TMP 
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APPENDIX B. TRANSFORMING THE CATEGORICAL RESPONSES TO 

AN INTERVAL SCALE 



A. RAW FREQUENCIES 



STATE VARIABLES RESPONSE FREQUENCIES 



QUES 

NUM 


RANGE 
(.METERS ) 


LOS 

(QTY ) 


SPEED 

(KMPH ) 


NO 

TMP 


MARG 

TMP 


EFF 

TMP 


GOOD 

TMP 


EXC 

TMP 


01 


3 


2400 


2 


2 . 5 


3 


22 . 0 


0 


9 


20 


11 


0 


02 


2 


1600 


4 


0.0 


i 


5.0 


0 


12 


19 


9 


0 


03 


4 


3200 


2 


2 . 5 


4 


40 . 0 


0 


2 


6 


20 


12 


04 


i 


750 


2 


2 . 5 


i 


5 . 0 


17 


18 


5 


0 


0 


05 


3 


2500 


3 


1 . 0 


i 


5 . 0 


3 


13 


20 


4 


0 


06 


2 


1400 


1 


4. 5 


2 


11 . 0 


5 


30 


5 


0 


0 


07 


4 


3400 


3 


1 . 0 


3 


24 . 0 


0 


0 


4 


22 


14 


08 


i 


640 


2 


2 . 5 


3 


24 . 0 


3 


22 


12 


3 


0 


09 


3 


2500 


1 


4. 5 


1 


5 . 0 


3 


27 


9 


1 


0 


10 


2 


1500 


3 


1 . 0 


3 


25.0 


0 


5 


20 


14 


1 


11 


4 


3300 


2 


2 . 5 


2 


12 . 0 


0 


9 


17 


12 


2 


12 


i 


800 


4 


0.0 


1 


5 . 0 


2 


13 


22 


3 


0 


13 


3 


2350 


3 


1 . 0 


4 


40 . 0 


0 


1 


8 


19 


12 


14 


2 


1575 


2 


2.5 


1 


5 . 0 


3 


30 


7 


0 


0 


15 


4 


3450 


3 


1 . 0 


4 


35.0 


0 


0 


5 


19 


16 


16 


1 


700 


1 


4 . 5 


2 


11 . 0 


16 


18 


6 


0 


0 


17 


3 


2450 


3 


1 . 0 


2 


10 . 0 


0 


3 


30 


7 


0 


18 


2 


1450 


2 


2 . 5 


4 


38.0 


0 


11 


23 


5 


1 


19 


4 


3400 


1 


4 . 5 


1 


5.0 


5 


20 


13 


1 


1 


20 


1 


575 


4 


0 . 0 


4 


42 . 0 


0 


3 


9 


12 


16 


21 


3 


2600 


4 


0.0 


3 


25.0 


0 


0 


6 


19 


15 


22 


2 


1525 


2 


2 . 5 


2 


11.0 


2 


22 


16 


0 


0 


23 


4 


3600 


4 


0 . 0 


1 


5 . 0 


2 


6 


16 


10 


6 


24 


1 


550 


2 


2 . 5 


4 


42.0 


2 


14 


16 


7 


1 


25 


3 


2625 


1 


4.5 


3 


22.0 


2 


19 


16 


3 


0 


26 


2 


1480 


4 


0 . 0 


4 


39.0 


0 


1 


8 


16 


15 


27 


4 


3300 


1 


4.5 


3 


25 . 0 


0 


10 


24 


5 


1 


28 


1 


660 


3 


1 . 0 


2 


10.0 


5 


18 


14 


3 


0 


29 


3 


2600 


4 


0 . 0 


1 


5 . 0 


0 


13 


19 


8 


0 


30 


2 


1585 


1 


4 . 5 


3 


23.0 


2 


30 


8 


0 


0 


31 


4 


3550 


4 


0 . 0 


3 


26.0 


0 


0 


3 


20 


17 


32 


1 


400 


1 


4 . 5 


3 


23.0 


14 


17 


7 


2 


0 


33 


3 


2640 


2 


2 . 5 


2 


11.0 


0 


19 


19 


2 


0 


34 


2 


1475 


3 


1 . 0 


1 


5.0 


2 


20 


18 


0 


0 


35 


4 


3740 


1 


4 . 5 


4 


45.0 


1 


5 


15 


17 


2 


36 


1 


760 


4 


0 . 0 


3 


20 . 0 


0 


3 


12 


22 


3 


37 


3 


2575 


4 


0.0 


4 


42 . 0 


0 


0 


0 


10 


30 


38 


2 


1480 


1 


4 . 5 


1 


5.0 


16 


22 


2 


0 


0 


39 


4 


3550 


3 


1.0 


1 


5.0 


1 


16 


13 


9 


1 


40 


1 


600 


1 


4.5 


1 


5.0 


25 


15 


0 


0 


0 


41 


3 


2305 


1 


4.5 


4 


40.0 


0 


15 


19 


6 


0 


42 


2 


1620 


4 


0.0 


2 


12.0 


0 


2 


23 


15 


0 
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43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 



1 

5 

0 

2 

1 

3 

1 

0 

2 

0 

0 

5 

1 

0 

1 

0 

35 

0 

0 

0 

9 

1 



4 


3300 


3 


1 . 0 


2 


10 . 0 


0 


8 


22 


9 


1 


575 


3 


1 . 0 


4 


42 . 0 


0 


7 


15 


13 


3 


2375 


2 


2.5 


1 


5.0 


5 


27 


8 


0 


2 


1565 


3 


1 . 0 


4 


40.0 


0 


6 


17 


15 


4 


3500 


1 


4.5 


2 


12 . 0 


3 


14 


20 


2 


1 


500 


4 


0.0 


2 


10 . 0 


0 


8 


15 


14 


3 


2480 


4 


0 . 0 


2 


11 . 0 


0 


4 


13 


22 


2 


1500 


1 


4.5 


4 


40 . 0 


0 


25 


13 


2 


4 


3450 


2 


2.5 


3 


21.0 


0 


10 


19 


9 


1 


700 


2 


2.5 


2 


10.0 


12 


22 


6 


0 


3 


2600 


3 


1 . 0 


3 


22.0 


0 


5 


22 


13 


2 


1390 


4 


0 . 0 


3 


24.0 


0 


0 


7 


28 


4 


3475 


2 


2.5 


1 


5 . 0 


5 


19 


13 


2 


1 


500 


1 


4.5 


4 


40 . 0 


10 


15 


10 


5 


3 


2475 


2 


2.5 


4 


42 . 0 


0 


9 


20 


10 


2 


1550 


3 


1 . 0 


2 


11.0 


0 


12 


25 


3 


4 


3600 


4 


0.0 


4 


40.0 


0 


0 


0 


5 


1 


540 


3 


1 . 0 


1 


5 . 0 


9 


20 


11 


0 


3 


2500 


1 


4.5 


2 


10 . 0 


9 


27 


4 


0 


2 


1500 


2 


2 . 5 


3 


24 . 0 


2 


18 


19 


1 


4 


3650 


4 


0 . 0 


2 


11 . 0 


0 


2 


13 


16 


1 


625 


3 


1.0 


3 


22.0 


3 


11 


17 


8 
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B. GROUPED FREQUENCIES 



STATE VARIABLES RESPONSE FREQUENCIES 



QUES 

NUH 


RANGE 
( METERS ) 


LOS 

(QTY ) 


SPEED 

( KMPH ) 


NO 

IMP 


MARG 

TUP 


EFF 

IMP 


GOOD 

IMP 


EXC 

IMP 


59 


4 


3600 


4 


0 . 0 


4 


40 . 0 


0 


0 


0 


5 


35 


37 


3 


2575 


4 


0 . 0 


4 


42.0 


0 


0 


0 


10 


30 


31 


4 


3550 


4 


0 . 0 


3 


26 . 0 


0 


0 


3 


20 


17 


07 


4 


3400 


3 


1 . 0 


3 


24.0 


0 


0 


4 


22 


14 


15 


4 


3450 


3 


1 . 0 


4 


35.0 


0 


0 


5 


19 


16 


21 


3 


2600 


4 


0 . 0 


3 


25.0 


0 


0 


6 


19 


15 


54 


2 


1390 


4 


0 . 0 


3 


24.0 


0 


0 


7 


28 


5 


13 


3 


2350 


3 


1 . 0 


4 


40.0 


0 


0 


9 


19 


12 


26 


2 


1480 


4 


0 . 0 


4 


39 . 0 


0 


0 


9 


16 


15 


20 


1 


575 


4 


0 . 0 


4 


42 . 0 


0 


3 


9 


12 


16 


03 


4 


3200 


2 


2 . 5 


4 


40.0 


0 


2 


6 


20 


12 


36 


i 


760 


4 


0 . 0 


3 


20.0 


0 


3 


12 


22 


3 


63 


4 


3650 


4 


0 . 0 


2 


11.0 


0 


2 


13 


16 


9 


44 


i 


575 


3 


1 . 0 


4 


42 . 0 


0 


7 


15 


13 


5 


11 


4 


3300 


2 


2.5 


2 


12 . 0 


0 


9 


17 


12 


2 


46 


2 


1565 


3 


1 . 0 


4 


40 . 0 


0 


6 


17 


15 


2 


48 


1 


500 


4 


0 . 0 


2 


10 . 0 


0 


8 


15 


14 


3 


51 


4 


3450 


2 


2 . 5 


3 


21 . 0 


0 


10 


19 


9 


2 


35 


4 


3740 


1 


4.5 


4 


45 . 0 


0 


6 


15 


17 


2 


23 


4 


3600 


4 


0 . 0 


i 


5 . 0 


0 


8 


16 


10 


6 


10 


2 


1500 


3 


1 . 0 


3 


25.0 


0 


5 


20 


15 


0 


43 


4 


3300 


3 


1 . 0 


2 


10 . 0 


0 


8 


22 


10 


0 


18 


2 


1450 


2 


2 . 5 


4 


38 . 0 


0 


11 


23 


6 


0 


27 


4 


3300 


1 


4.5 


3 


25.0 


0 


10 


24 


6 


0 


49 


3 


2480 


4 


0 . 0 


2 


11.0 


0 


4 


13 


23 


0 


57 


3 


2475 


2 


2.5 


4 


42 . 0 


0 


9 


20 


11 


0 


39 


4 


3 550 


3 


1 . 0 


1 


5 . 0 


0 


17 


13 


10 


0 


42 


2 


1620 


4 


0 . 0 


2 


12 . C 


0 


2 


23 


15 


0 


53 


3 


2600 


3 


1 . 0 


3 


22 . 0 


0 


5 


22 


13 


0 


01 


3 


2400 


2 


2 . 5 


3 


22 . 0 


0 


9 


20 


11 


0 


29 


3 


2600 


4 


0.0 


1 


5 . 0 


0 


13 


19 


8 


0 


02 


2 


1600 


4 


0.0 


1 


5.0 


0 


12 


19 


9 


0 


17 


3 


2450 


3 


1.0 


2 


10.0 


0 


3 


30 


7 


0 


41 


3 


2305 


1 


4.5 


4 


40.0 


0 


15 


19 


6 


0 


50 


2 


1500 


1 


4.5 


4 


40 . 0 


0 


25 


13 


2 


0 


58 


2 


1550 


3 


1 . 0 


2 


11.0 


0 


12 


25 


3 


0 


33 


3 


2640 


2 


2 . 5 


2 


11.0 


0 


19 


19 


2 


0 


24 


1 


550 


2 


2 . 5 


4 


42 . 0 


2 


14 


16 


8 


0 


64 


1 


625 


3 


1 . 0 


3 


22.0 


3 


11 


17 


9 


0 


47 


4 


3500 


1 


4.5 


2 


12.0 


3 


14 


20 


3 


0 


55 


4 


3475 


2 


2 . 5 


1 


5 . 0 


5 


19 


13 


3 


0 


19 


4 


3400 


1 


4.5 


1 


5 . 0 


5 


20 


13 


2 


0 


05 


3 


2500 


3 


1 . 0 


1 


5 . 0 


3 


13 


20 


4 


0 


56 


1 


500 


1 


4.5 


4 


40 . 0 


10 


15 


10 


5 


0 


32 


1 


400 


1 


4.5 


3 


23.0 


14 


17 


7 


2 


0 


25 


3 


2625 


1 


4.5 


3 


22 . 0 


2 


19 


16 


3 


0 


08 


1 


640 


2 


2 . 5 


3 


24 . 0 


3 


22 


12 


3 


0 


12 


1 


800 


4 


0 . 0 


1 


5.0 


2 


13 


22 


3 


0 


28 


1 


660 


3 


1 . 0 


2 


10.0 


5 


18 


14 


3 


0 
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62 

09 

34 

22 

14 

60 

52 

45 

30 

16 

06 

04 

61 

38 

40 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



2 


1500 


2 


2 . 5 


3 


24.0 


2 


18 


3 


2500 


1 


4.5 


1 


5 . 0 


3 


27 


2 


1475 


3 


1 . 0 


1 


5 . 0 


2 


20 


2 


1525 


2 


2 . 5 


2 


11 . 0 


2 


22 


2 


1575 


2 


2 . 5 


1 


5 . 0 


3 


30 


1 


540 


3 


1.0 


1 


5.0 


9 


20 


1 


700 


2 


2.5 


2 


10 . 0 


12 


22 


3 


2375 


2 


2 . 5 


1 


5 . 0 


5 


27 


2 


1585 


1 


4.5 


3 


23 . 0 


2 


30 


1 


700 


1 


4.5 


2 


11.0 


16 


18 


2 


1400 


1 


4.5 


2 


11 . 0 


5 


30 


1 


750 


2 


2.5 


1 


5 . 0 


17 


18 


3 


2500 


1 


4.5 


2 


10 . 0 


9 


27 


2 


1480 


1 


4.5 


1 


5 . 0 


16 


22 


1 


600 


1 


4.5 


1 


5 . 0 


25 


15 



20 

10 

18 

16 

7 

11 



53 
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C. RELATIVE FREQUENCIES 



RESPONSE FREQUENCIES 



EXC 

IMP 


GOOD 

IMP 


EFF 

TUP 


MARC 

IMP 


NO 

TMP 


0.875 


0.125 


0 


0 


0 


0.75 


0.25 


0 


0 


0 


0.425 


0 . 5 


0.075 


0 


0 


0.35 


0 . 55 


0.1 


0 


0 


0 . 4 


0.475 


0.125 


0 


0 


0.375 


0.475 


0.15 


0 


0 


0.125 


0 . 7 


0.175 


0 


0 


0 . 3 


0.475 


0.225 


0 


0 


0.375 


0.4 


0.225 


0 


0 


0 • 4 


0.3 


0.225 


0 . 075 


0 


0 . 3 


0. 5 


0.15 


0.05 


0 


0.075 


0. 55 


0 . 3 


0.075 


0 


0.225 


0.4 


0.325 


0.05 


0 


0 .125 


0.325 


0 .375 


0 .175 


0 


0.05 


0.3 


0.425 


0.225 


0 


0 .05 


0.375 


0 .425 


0.15 


0 


0 . 075 


0.35 


0 .375 


0 . 2 


0 


0.05 


0.225 


0.475 


0.25 


0 


0 .05 


0.425 


0 .375 


0.15 


0 


0 .15 


0 .25 


0.4 


0 . 2 


0 


0 


0.375 


0.5 


0.125 


0 


0 


0.25 


0.55 


0.2 


0 


0 


0.15 


0.575 


0.275 


0 


0 


0.15 


0.6 


0 . 25 


0 


0 


0.575 


0.325 


0.1 


0 


0 


0.275 


0 . 5 


0.225 


0 


0 


0.25 


0.325 


0.425 


0 


0 


0.375 


0.575 


0.05 


0 


0 


0.325 


0.55 


0.125 


0 


0 


0.275 


0 . 5 


0 .225 


0 


0 


0 . 2 


0.475 


0.325 


0 


0 


0.225 


0.475 


0 . 3 


0 


0 


0.175 


0.75 


0.075 


0 


0 


0 .15 


0 .475 


0.375 


0 


0 


0.05 


0.325 


0.625 


0 


0 


0 . 075 


0 .625 


0 . 3 


0 


0 


0.05 


0.475 


0.475 


0 


0 


0.2 


0 . 4 


0.35 


0.05 


0 


0.225 


0.425 


0.275 


0.075 


0 


0.075 


0. 5 


0.35 


0.075 


0 


0 . 075 


0.325 


0.475 


0 .125 


0 


0.05 


0.325 


0 . 5 


0 .125 


0 


0 . 1 


0 . 5 


0.325 


0.075 


0 


0.125 


0 .25 


0 .375 


0.25 


0 


0.05 


0.175 


0.425 


0.35 


0 


0.075 


0 . 4 


0 .475 


0 .05 


0 


0 . 075 


0 . 3 


0 . 55 


0.075 


0 


0.075 


0 .55 


0.325 


0.05 


0 


0 . 075 


0.35 


0.45 


0.125 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0.5 


0.45 


0 .05 


0.25 


0.675 


0.075 


0.45 


0.5 


0.05 


0.4 


0.55 


0.05 


0.175 


0.75 


0.075 


0.275 


0.5 


0.225 


0.15 


0.55 


0 . 3 


0.2 


0.675 


0.125 


0.2 


0 .75 


0.05 


0.15 


0.45 


0.4 


0.125 


0.75 


0.125 


0.125 


0.45 


0.425 


0.1 


0.675 


0.225 


0 . 05 


0.55 


0.4 


0 


0.375 


0.625 
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D. CUMULATIVE RELATIVE FREQUENCIES 



RESPONSE FREQUENCIES 



EXC 

IMP 


GOOD 

IMP 


EFF 

IMP 


MARG 

TMP 


NO 

TMP 


0.875 


1 


1 


1 


1 


0.75 


1 


1 


1 


1 


0.425 


0.925 


1 


1 


1 


0.35 


0 . 9 


1 


1 


1 


0.4 


0.875 


1 


1 


1 


0.375 


0.85 


1 


1 


1 


0.125 


0.825 


1 


1 


1 


0 . 3 


0.775 


1 


1 


1 


0 .375 


0.775 


1 


1 


1 


0 • 4 


0 . 7 


0.925 


1 


1 


0.3 


0 . 8 


0 .95 


1 


1 


0.075 


0.625 


0.925 


1 


1 


0.225 


0.625 


0 .95 


1 


1 


0.125 


0 .45 


0 .825 


1 


1 


0.05 


0.35 


0 . 775 


1 


1 


0.05 


0.425 


0.85 


1 


1 


0 . 075 


0.425 


0 . 8 


1 


1 


0.05 


0 . 275 


0 .75 


1 


1 


0.05 


0.475 


0.85 


1 


1 


0.15 


0 . 4 


0 . 8 


1 


1 


0 


0.375 


0.875 


1 


1 


0 


0.25 


0 . 8 


1 


1 


0 


0.15 


0.725 


1 


1 


0 


0 .15 


0.75 


1 


1 


0 


0 . 575 


0.9 


1 


1 


0 


0.275 


0.775 


1 


1 


0 


0.25 


0 . 575 


1 


1 


0 


0 .375 


0.95 


1 


1 


0 


0.325 


0 . 875 


1 


1 


0 


0.275 


0.775 


1 


1 


0 


0.2 


0.675 


1 


1 


0 


0 . 225 


0.7 


1 


1 


0 


0 .175 


0.925 


1 


1 


0 


0.15 


0.625 


1 


1 


0 


0 .05 


0.375 


1 


1 


0 


0 .075 


0.7 


1 


1 


0 


0.05 


0.525 


1 


1 


0 


0 . 2 


0.6 


0 .95 


1 


0 


0.225 


0.65 


0.925 


1 


0 


0.075 


0 . 575 


0.925 


1 


0 


0.075 


0.4 


0.875 


1 


0 


0.05 


0.375 


0.875 


1 


0 


0.1 


0.6 


0.925 


1 


0 


0 . 125 


0.375 


0.75 


1 


0 


0.05 


0.225 


0.65 


1 


0 


0.075 


0.475 


0.95 


1 


0 


0 .075 


0.375 


0.925 


1 


0 


0 .075 


0.625 


0.95 


1 


0 


0 .075 


0.425 


0.875 


1 
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0 


0 


0 . 5 


0 .95 


1 


0 


0 


0.25 


0 .925 


1 


0 


0 


0.45 


0 .95 


1 


0 


c 


0.4 


0 .95 


1 


0 


0 


0.175 


0 . 925 


1 


0 


0 


0.275 


0.775 


1 


0 


0 


0 .15 


0 . 7 


1 


0 


0 


0 . 2 


0.875 


1 


0 


0 


0 . 2 


0.95 


1 


0 


0 


0 .15 


0 . 6 


1 


0 


0 


0.125 


0.875 


1 


0 


0 


0.125 


0 . 575 


1 


0 


0 


0 . 1 


0.775 


1 


0 


0 


0.05 


0.6 


1 


0 


0 


0 


0.375 


1 



E. EXTREME INSTANCES 
EXCELLENT CASES 

~C7S75 - I I 1 I 

0.75 1111 

POOR CASES 

__ _ _ _______ 
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F. GROUPED CUMULATIVE RELATIVE FREQUENCY MATRICES 

CEO UP I 



0.425 


0 .925 


1 


1 


1 


0.35 


0 . 9 


1 


1 


1 


0.4 


0.875 


1 


1 


1 


0.375 


0.85 


1 


1 


1 


0.125 


0 .825 


1 


1 


1 


0 . 3 


0.775 


1 


1 


1 


0.375 


0.775 


1 


1 


1 


TGROUPI 






0.425 


0.925 








0.35 


0 . 9 








0 . 4 


0.875 








0.375 


0.85 








0 .125 


0.825 








0 . 3 


0.775 








0.375 


0.775 








GROUPII 






0.4 


0 . 7 


0.925 


1 


1 


0 . 3 


0 . 8 


0.95 


1 


1 


0.075 


0.625 


0.925 


1 


1 


0 .225 


0.625 


0.95 


1 


1 


0.125 


0.45 


0.825 


1 


1 


0.05 


0.35 


0.775 


1 


1 


0.05 


0.425 


0.85 


1 


1 


0.075 


0.425 


0 . 8 


1 


1 


0.05 


0.275 


0.75 


1 


1 


0.05 


0.475 


0.85 


1 


1 


0.15 


0.4 


0 . 8 


1 


1 


TGROUPII 






0.4 


0.7 


0.925 






0 . 3 


0 . 8 


0.95 






0 .075 


0.625 


0.925 






0.225 


0.625 


0.95 






0.125 


0.45 


0.825 






0.05 


0.35 


0.775 






0.05 


0.425 


0.85 






0 .075 


0.425 


0 . 8 






0 .05 


0.275 


0.75 






0 .05 


0.475 


0.85 






0.15 


0 . 4 


0 . 8 
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GROUPIII 



0.375 


0.875 


1 


1 


0.25 


0 . 8 


1 


1 


0.15 


0 .725 


1 


1 


0.15 


0.75 


1 


1 


0.575 


0.9 


1 


1 


0.275 


0.775 


1 


1 


0.25 


0.575 


1 


1 


0.375 


0.95 


1 


1 


0.325 


0.875 


1 


1 


0.275 


0.775 


1 


1 


0.2 


0.675 


1 


1 


0 .225 


0.7 


1 


1 


0.175 


0.925 


1 


1 


0 . 15 


0.625 


1 


1 


0 .05 


0.375 


1 


1 


0.075 


0.7 


1 


1 


0.05 


0.525 


1 


1 



TGROUPIII 



0.375 


0.875 


0.25 


0 . 8 


0.15 


0.725 


0.15 


0.75 


0.575 


0 . 9 


0.275 


0.775 


0.25 


0.575 


0.375 


0.95 


0 .325 


0 .875 


0.275 


0.775 


0.2 


0.675 


0 .225 


0.7 


0.175 


0.925 


0.15 


0.625 


0 .05 


0.375 


0.075 


0.7 


0.05 


0.525 
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oooooooooooo 



CROUPIV 



0 . 2 


0 . 6 


0 .95 


1 


0 .225 


0 .65 


0 . 925 


1 


0 .075 


0 .575 


0 . 925 


1 


0 .075 


0.4 


0.875 


1 


0 .05 


0.375 


0 .875 


1 


0 . 1 


0 . 6 


0 .925 


1 


0.125 


0.375 


0.75 


1 


0.05 


0 . 225 


0.65 


1 


0.075 


0.475 


0 .95 


1 


0.075 


0.375 


0 . 925 


1 


0.075 


0.625 


0.95 


1 


0 .075 


0.425 


0.875 


1 



TCROUPIV 



0 . 2 


0.6 


0.95 


0 .225 


0.65 


0 .925 


0.075 


0.575 


0 .925 


0.075 


0.4 


0.875 


0.05 


0.375 


0 .875 


0 . 1 


0 . 6 


0.925 


0.125 


0.375 


0.75 


0.05 


0 . 225 


0 .65 


0.075 


0.475 


0 .95 


0.075 


0.375 


0 . 925 


0.075 


0.625 


0 .95 


0.075 


0.425 


0 .875 
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GROUPV 



0 0.5 0.95 1 
0 0.25 0.925 1 
0 0.45 0.95 1 
0 0.4 0.95 1 
0 0.175 0.925 1 
0 0.275 0.775 1 
0 0.15 0.7 1 
0 0.2 0.875 1 
0 0.2 0.95 1 
0 0.15 0.6 1 
0 0.125 0.875 1 
0 0 . 1250.5751 
0 0.1 0.7751 
0 0.05 0.6 1 

T GROUPV 



0.5 


0.95 


0 .25 


0.925 


0.45 


0.95 


0.4 


0 .95 


0.175 


0 .925 


0.275 


0.775 


0.15 


0.7 


0.2 


0.875 


0.2 


0 .95 


0.15 


0 . 6 


0.125 


0.875 


0.125 


0.575 


0.1 


0.775 


0.05 


0.6 
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G. COM PITER PROGRAMS USED 

THE APL COMPUTER PROGRAM NORM 



V NORM l DDV 

V NORM 

[ 1] *VNORM\MTX\NOR\ROWAV\GRAV\S\AA\B\AAI\II 

[2] nV NORM ;MTX ; NOR; ROW AV-,GRAV;S;AA;B;AAI; II 

[3] □-«- 1 INPUT THE CUMREL FREQUENCY MATRIX' 

[4] MTX + □ 

[5] NOR+NQUAN MTX 

[6] S+pMTX 

[7] R0WAV+(+/N0R')*(SZ2l) 

[8] ROWAV+iSZll ,l)pROWAV 

[9] COLAV+ ( +/NOR ) t (S [ 1 ] ) 

[10] GRAV+( + / ( + /NOR))i USZ 1] )x(S[2] )) 

[11] □ +' NORMALIZED VALUES ROW AVERAGE' 

[ 12 ] □<-' ' 

[13] □ +NOR,ROWAV 

[14] □«-' ' 

[15] ' COLUMN AVERAGES ' 

[16] ' 

[17] U+COLAV 

[18] □<-« ' 

[19] □ +' GRAND AVERAGE' 

[ 20 ] 1 ' 

[21] U+GRAV 

[22] Ck 1 ' 

[23] AAI+SpQ 

[24] 11 + 0 

[25] L2 '.II+II+l 

[26] AA+(NORh (ID!- ,ROWAV )* 2 

[27] AAIL;II1+AA 

[28] -KIKS[2] )/L2 

[29] AAI+(S)pAAI 

[30] AAI++/AAI 

[31] B+ + / (.(COLAV-GRAV)* 2) 

[32] SQR+(BiAAI)* 0.5 

[33] SQR+iiSZll ) ,1 )pSQR 

[34] SSI+GRAV-(RONAVxSQR ) 

[35] ^(Stl],l)pMI 

[36] GRAV+CSZll t l)pGRAV 

[37] B ' 

[38] D+ ' ' 

[39] □ +B 

[40] ' ' 

[41] ' AI ' 

[42] ' 

[43] Q+AAI 

[44] ' 

[45] □ +' SCALE VALUES = GRAND AVERAGE - ( ROW AVERAGE x (B+AI) * 

[46] □-«- ' 

[47] U+SSI ,GRAV ,ROWAV ,SQR 

[48] ' ' 

[49] D+ ' COLUMN AVERAGES ' 

[50] 0+ ' ' 

[51] 0 +COLAV 

V 
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VNQUAN CD]V 

V Z*NQUAN P;A;B',C',D 

[1] *( (+/U+(P<0 )v(p*i ) ) )>o )/Ll 

[2] C«- 2.515517 0.802853 0.010328 

[3] D* 1.432788 0.189269 0.001308 

[4] P*UA*(P< 0 . 5 ) )xp) + ( (p>0 . 5 )x (l-P) ) 

C 5] B*(®P* 2)*0.5 

[6] Z«-((2xA)-i)x-E-((Bo,* o 1 2 )+ . xC)* (l + ( (B® . * l 2 3)+.*L)) 

[7] +0 

[8] Ll:U*' THERE IS NO QUANTILE FOR P - ' ,$A/P 

V 



THE APL COMPUTER PROGRAM TRANS 



V TRANS CD] 7 

V TRANS 

[ 1 ] o V TRANS ; UP ; LOW ; BETA ; ALPHA ; MX ; COL UP ; COLLOW ; TRCOL 

[2] !> ' IJVPC/Z’ COL AVERAGES' 

[3] CCL^O 

[4] S^-pCCL 

C 5] COLUP*COL [S] 

[6] COLLOW*COL [1] 

[7] □«- ' IJVPi/Z 1 THE VECTOR OF BOUNDRY VALUES TO BE TRANSFORMED ' 

[8] MX * □ 

[9] D*' INPUT UPPER BOUNDRY VALUE' 

[10] UP * □ 

[11] D*' INPUT LOWER BOUNDRY VALUE' 

[12] LOW * □ 

[13] BETA.* ( UP -LOW ) * (CO LUP- COLLOW ) 

[14] A£PffA«-0P-flEMxC0L0P 

[15] TR*ALPHA+BETA*MX 

[16] TRCOL*ALPHA+BETA*COL 

[17] 1 TRANSFORMED UPPER BOUNDS ARE' 

[18] OPPCCL 

[19] ' 

[20] 0-e- * TRANSFORMED DATA ' 

[21] □ *TR 

V 



63 



H. THE NORMALIZED AND TRANSFORMED VALUES OF THE INSTANCES 



T GROUP I I 

NORMALIZED VALUES RON AVERAGE 



"0 .2529332678 
"0 .5240018704 

1 .43980047 
0 .7551784916 
“1 .150435626 

1 .64521144 
“1 . 64521144 
"1.43980047 
"1 .64521144 

1 .64521144 
"1 . 036431485 



0 . 5240018704 
0 . 8414567174 
0.3181998625 
0.3181998625 
0.1253809931 
0.384877085 
0 .1887560404 
0 .1887560404 
0 . 5974048985 
0 . 06254483635 
0 .2529332678 



1.43980047 

1 . 64521144 

1.43980047 

1.64521144 
0.9345033954 
0 .7551784916 
1 . 036431485 
0.8414567174 
0 .67418914 
1 . 036431485 
0 . 8414567174 



0.5702896907 
0.6542220957 
0 .1060666208 
_0 . 4027442703 
"0 .1137710747 
_0 .4249700112 
"0 . 2658453318 
"0 .2623665976 
"0 . 5228090662 
_0. 2237749304 
”0 . 1493026786 



COLUMN AVERAGES 

1.198129767 0.01829137737 1.117242841 



GRAND AVERAGE 
"0 .02086518298 



B 



2 .682775012 



AI 



1.435887596 
2 .405328399 
4 .213551553 
2.891657548 
2 .173687547 
2 .883347076 
3.604518528 
2 .610194948 
2.698156333 
3 .634597104 
1 .779340999 

SCALE VALUES = 



0 . 8003854753 
”0.7117889701 
0 .1054995408 
0.4087904496 
0 . 1055285172 
0 .3890574987 
0.2084841776 
0 . 2451241358 
0 . 5004515703 
0.1713888323 
0.1624633793 



GRAND AVERAGE - 



”0.02086518298 
0. 02086518298 
”0 .02086518298 
_0 . 02086518298 
0 . 02086518298 
0.02086518298 
0.02086518298 
0. 02086518298 
”0 .02086518298 
0.02086518298 
0.02086518298 



( ROW AVERAGE x 



0.5702896907 
0 .6542220957 
0 .1060666208 
0.4027442703 
"0.1137710747 
0.4249700112 
0.2658453318 
”0.2623665976 
_0 . 5228090662 
0.2237749304 
0.1493026786 



(B*AI) * .5) 



1.366884769 
1 .056099743 
0.7979358365 
0 . 9632049299 
1 .110947581 
0 . 9645920205 
0 . 8627172766 
1.01380786 
0 .9971455871 
0 .8591400964 
1.227898683 



64 



COLUMN AVERAGES 



"1.198129767 0.01829137737 1.117242841 

TRANS 

INPUT UPPER BOUNDRY VALUE 

□ : 

100 . 0 

INPUT LOWER BOUNDRY VALUE 

□ : 

0 . 0 

TRANSFORMED UPPER BOUNDS ARE 
1 . 105427358E“15 53.53673384 100 

TRANSFORMED DATA 
17.17841399 
21 . 00486097 
47 .1902545 
34.09124367 
56 .30447038 
68 . 54997163 
60 .75108343 
62 .3335483 
73 . 36103621 
59 . 14895057 
58 . 76346389 



T GROUP I 

NORMALIZED VALUES ROW AVERAGE 



0 .1887560404 
~ 0 .384877085 
0.2529332678 
“0 .3181998625 
“1.150435626 
_0 .5240018704 
“0 .3181998625 



1.43980047 
1 .281728757 
1 .150435626 
1 .036431485 
0 .9345033954 
0 .7551784916 
0 .7551784916 



0.6255222146 
0 . 4484258358 
0.4487511792 
0 .3591158114 
“0 .1079661155 
0 .1155883106 
0 . 2184893146 



COLUMN AVERAGES 
“0.4482005164 1.050465245 
GRAND AVERAGE 
0.3011323644 



B 



1 .122999532 



AI 



1 .326098153 
1 .388787515 
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0 . 9847221265 
0 .917513044 
2 . 173485362 
0 .8181511992 
0 . 5760705455 

SCALE VALUES = GRAND AVERAGE - ( ROW AVERAGE x (BiAI) * .5) 



0.2744993195 

0.1021066199 

0 .1780916163 

0.09616729347 

0.3787390085 

0 .1657110832 

0.003925450303 



0.3011323644 
0 .3011323644 
0 .3011323644 
0.3011323644 
0 .3011323644 
0 .3011323644 
0 .3011323644 



0.6255222146 
0.4484258358 
0.4487511792 
0 .3591158114 
0 .1079661155 
0.1155883106 
0 .2184893146 



0.9202417922 
0 .8992322745 
1 . 067905786 
1.106327389 
0 .7188055603 
1 .17158284 
1.396213885 



COLUMN AVERAGES 



0.4482005164 1.050465245 
TRANS 

INPUT UPPER BOUNDRY VALUE 

□ : 

53 . 536734 

INPUT LOWER BOUNDRY VALUE 

□ : 

0 . 0 

TRANSFORMED UPPER BOUNDS ARE 
3 . 5527136 7 9E 15 53.536734 



TRANSFORMED DATA 
6.20511592 
12.36348848 
9 . 649081675 
12.57565863 
29.54070381 
21.93072188 
15.87080765 



TGROUPIII 

NORMALIZED VALUES ROW AVERAGE 



"0 .3181998625 
“0.67418914 
1.036431485 
1 .036431485 
0.1887560404 
0 .5974048985 
0.67418914 
0.3181998625 
0.4533333679 
0 . 5974048985 
0 . 8414567174 
0.7551784916 
“0 . 9345033954 



1 .150435626 
0 . 8414567174 
0 . 5974048985 
0.67418914 
1.281728757 
0.7551784916 
0 .1887560404 
1.64521144 
1.150435626 
0 .7551784916 
0.4533333679 
0 . 5240018704 
1.43980047 



0.4161178819 
0.08363378866 
0 .2195132933 
“0 . 1811211726 
0 . 7352423985 
0. 07888679654 
0 .2427165498 
0.6635057888 
0.3485511292 
0. 07888679654 
0 .1940616747 
0.1155883106 
0 .2526485371 
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1 . 036431485 
1 .64521144 
1.43980047 
1 . 64521144 



0.3181998625 0.3591158114 
0.3181998625 "0.9817056513 
0.5240018704 0.4578992996 
0.06254483635 0.7913333019 



COLUMN AVERAGES 
0.8126365611 0.7084504496 



GRAND AVERAGE 
"0.05209305576 



B 



1 .156852847 



AI 



1 .078445099 
1 . 148591183 
1 .334710664 
1 . 463111462 
0 . 5972946791 
0 .9147409136 
0 .3723371922 
1 . 927491972 
1 . 286037493 
0.9147409136 
0 .8382406825 
0 .8181511992 
2 .818659422 
0 . 91751304U 
0 . 8804798637 
1 .928259815 
1.45821575 

SCALE VALUES = 



_0. 4830723627 
0 .1360270889 
0.1522718504 
0 .1089602245 
_1 .075327804 
0 . 1408074984 
0.3757359646 
0 . 5661218971 
_0 .3826747026 
0 . 1408074984 
0 .1758857464 
0.08535423712 
0.2139511961 
0 .3511505256 
1 . 07318704 
0 .3025784214 
0 . 6 527421S52 



GRAND AVERAGE - 



0 .05209305576 
0.05209305576 
0 .05209305576 
0 .05209305576 
0 .05209305576 
"0.05209305576 
0.05209305576 
"0.05209305576 
0.05209305576 
”0.05209305576 
0.05209305576 
0.05209305576 
"0 .05209305576 
0. 05209305576 
_0 .05209305576 
0.05209305576 
0 .05209305576 



(ROW AVERAGE x 



0.4161178819 
0 .08363378866 
_0 .2195132933 
0 .1811211726 
0 .7352423985 
_0 .07888679654 
0.2427165498 
0.6635057888 
0.3485511292 
_0 .07888679654 
0 .1940616747 
0 .1155883106 
0 .2526485371 
"0.3591158114 
_0 .9817056513 
0.4578992996 
0 . 7913333019 



(B±AI) * .5) 



1 .035714459 
1 . 00358999 
0 .9309910257 
0 .8892018419 
1.391697147 
1 .124579099 
1 . 762669339 
0 .774716438 
0 . 9484452039 
1 .124579099 
1 .174774991 
1 .189110665 
0 .640645468 
1 .122878939 
1 . 146249992 
0.7745621745 
0.8906932632 
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COLUMN AVERAGES 



0.8126365611 0.7084504496 

TRANS 

INPUT UPPER BOUNDRY VALUE 

□ : 

100.0 

INPUT LONER BOUNDRY VALUE 

□ : 

53 . 536734 

TRANSFORMED UPPER BOUNDS ARE 
53.536734 100 

TRANSFORMED DATA 
63.60363281 
74 . 20451018 
83.01091652 
81.68791555 
45.51254274 
74.05848734 
89 . 83687224 
61.06679397 
66 .67038912 
74.05848734 
83 .7322288 
80 . 96684485 
71 .82423312 
89 .08588313 
111 . 141278 
87.60219509 
98.29833015 



BET*- (10 0-53 . 536734 )t(1.329835441-. 06 802605834 ) 

BET 

36 .82272983 

ALP<rl 0 0 + (BETx . 06802605834 ) 

ALP 

102 . 5049052 

UPE+ALP+iBET* 1 . 267100775) 

UPE 

149.1630147 



T GROUP IV 

NORMALIZED VALUES ROW AVERAGE 



0 . 8414567174 
"0 .7551784916 
"1 .43980047 
“1 .43980047 



0.2529332678 
0 .384877085 
0 .1887560404 
0 .2529332678 



1 .64521144 
1.43980047 
1.43980047 
1 .150435626 



0 .3522293302 
0.3564996877 
0 . 06291868014 
“0.1807660371 
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“1 .64521144 
1 . 281728757 
“1.150435626 
1 .64521144 
1 .43980047 
_1 .43980047 
1 . 43980047 
1 .43980047 



0 .3181998625 
0.2529332678 
0 .3181998625 
0.7551784916 
0 . 06254483635 
0.3181998625 
0.3181998625 
0 . 1887560404 



1 .150435626 
1.43980047 
0 .67418914 
0 .384877085 
1 .64521144 
1.43980047 
1 .64521144 
1 .150435626 



'0 . 2709918921 
_0 .1370016603 
0 . 2648154496 
"0.6718376156 
0 . 04762204472 
"0.1060666208 
0 . 1745369443 
"0 .1593736279 



COLUMN AVERAGES 

1.329835441 "0.06802605834 1.267100775 



GRAND AVERAGE 
"0 . 04358690797 



B 



3 .372934763 



AI 



3 .106548825 
2.410174235 
4 . 169803347 
3.36247368 
3 .911164149 
3 . 723520871 
1.668902612 
2.071048262 
4 .776854354 
4.213551553 
4 .789607793 
3.355956505 

SCALE VALUES = GRAND AVERAGE - ( ROW AVERAGE x (B*AI) * .5) 



_0 .4106074699 
0.4653211892 
0 . 100175096 
0 .1374601038 
0 .2080690671 
"0.1739794959 
0.3328843213 
0 .8137927532 
0.08360358382 
0.05131141717 
0 .1900544398 
0.1161893582 

COLUMN AVERAGES 



_0 . 043 58690797 
0 . 04358690797 
0 . 04358690797 
"0 . 04358690797 
"0 . 04358690797 
0 . 04358690797 
"0 . 04358690797 
0 . 04358690797 
"0.04358690797 
"0.04358690797 
"0.04358690797 
"0 . 04358690797 



0.3522293302 
0.3564996877 
_0 . 06291868014 
"0 .1807660371 
"0 . 2709918921 
0.1370016603 
"0.2648154496 
0.6718376156 
0 . 04762204472 
"0.1060666208 
0 .1745369443 
"0.1593736279 



1.041993186 
1 .182986398 
0 .8993861263 
1.001554356 
0 .9286476179 
0 .9517591802 
1.4216362 
1 .276170969 
0 . 8402973052 
0 . 8947048978 
0 .8391778165 
1 .002526379 



"1.329835441 0.06802605834 1.267100775 

TRANS 
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INPUT UPPER BOUNDRY VALUE 

□ : 

UPE 

INPUT LOWER BOUNDRY VALUE 

□ : 

53 . 536734 

TRANSFORMED UPPER BOUNDS ARE 
53.536734 99.99999999 149.1630147 

TRANSFORMED DATA 
87.38521723 
85.37050872 
98.81618467 
107 . 5665614 
110 . 1665762 
96.09850519 
114.7626146 
132.4709759 
99.42639298 
104 .3943316 
95 . 50658187 
106 . 7833145 



TGROUPV 

NORMALIZED VALUES ROW AVERAGE 



1 . 010 0667 5 7E_7 
"6 . 741891400E_1 
1.25380993 1E~ 1 
"2 . 52933267 8E_ 1 
9 . 3 4 5 0 3 3 9 54E - 1 
“5 . 974048985E 1 
"1 . 036431485E0 
8 . 414567174E'l 
_8 . 414567 174E 1 
“1 . 03643148 5E0 
“1 . 150435626E0 
"1 . 15 043 5626E0 
'1 . 281728757E0 
"1 . 64521 144 0E0 



1 . 64 5 2 1 144 0E0 
1 .439800470E0 
1.645211440E0 
1.645211440E0 
1 . 439800470E0 
7 . 5 517 84916E_1 
5.240018704E 1 
1 . 150435626E0 
1 . 64521 144 0E0 
2 . 529332678E'l 
1 . 150435626E0 
1 . 887560404E“1 
7 . 551784916E_1 
2.5293 3 26 78E 1 



8 . 226057706E"1 
3.8280 5664 8£ 1 
7 . 599152235E'l 
6 . 9613 90 86 2E~ 1 
2 . 526485371E 1 
_7 . 888679654E_2 
2 . 562148074E 1 
1. 544894545E_1 
4 . 018773614E_1 
”3 . 917491087E 1 
0 . 000000000E0 
4.808397 92 9E 1 
”2.632751325E_1 
”6 . 961390862E 1 



COLUMN AVERAGES 
0.8262856749 1.035021384 



GRAND AVERAGE 
0 .1043678548 



B 

1 .732231984 
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AT 



1 .353360175 
2 .234476035 
1 . 567498782 
1 .801476666 
2 . 818659422 
0 .9147409136 
1 .217476129 
1 .983817554 
3 . 091759263 
0 .8312307332 
2 .64700426 
0 .8967171601 
2 .074495569 
1 . 801476666 

SCALE VALUES - GRAND AVERAGE - ( ROW AVERAGE * (B*AI ) 



0 . 8262857891 
0 .2326815247 
0 .6944809054 
0 . 5782611302 
0 .09369293795 
0 .004189260708 
0 .4099845144 
0 .0399935175 
0 .1964431591 
0 .6698908224 
0 .1043678548 
0.7726743588 
0 .3449462429 
0 .7869968397 



0 .1043678548 
0 .1043678548 
0 .1043678548 
0 .1043678548 
0 .1043678548 
0 .1043678548 
0 . 1043678548 
0 .1043678548 
0 .1043678548 
0 .1043678548 
0 . 1043678548 
0.1043678548 
0 .1043678548 
0 .1043678548 



0.8226057706 1 . 

0.3828056648 0 . 

0.7599152235 1 . 

0.6961390862 0 . 

0.2526485371 0 . 

_ 0 . 07888679654 1 . 

0.2562148074 1 . 

0.1544894545 0 . 

0.4018773614 0 . 

0.3917491087 1 . 

0 0 . 

0.4808397929 1 . 

_ 0 . 2632751325 0 . 

0.6961390862 0 . 



COLUMN AVERAGES 



" 0.8262856749 1.035021384 

TRANS 

INPUT UPPER BOUNDRY VALUE 

□ : 

UPE 

INPUT LONER BOUNDRY VALUE 

0 : 

100.0 

TRANSFORMED UPPER BOUNDS ARE 
100 149.1630147 



TRANSFORMED DATA 
99.99999698 
115.6789657 
103.4813814 
106 .5511138 
119.3500945 
121 .7141701 
132.6538113 
120 . 7684663 



* . 5 ) 



131348305 

880471243 

051234053 

9805928133 

7839380152 

37611261 

192814197 

9344415953 

7485144544 

443584567 

8089575013 

389873538 

9137907782 

9805928133 
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116.6361357 
139 . 5187601 
124.5815073 
142 .2335988 
130.9359445 
142.6119009 
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I. THE REGRESSION MATRIX 





STATE VARIABLES 






SCALE VALUE 


QUES 

NUM 


RANGE 

( METERS ) 


LOS X 
1000 (QTY) 


SPEED X 
100 ( KMPH ) 


TMP 

RATING 


31 


3550 


0 


2600 


6 . 20511592 


7 


3400 


1000 


2400 


12 . 36348848 


15 


3450 


1000 


3500 


9 . 64908167 


21 


2600 


0 


2500 


12 . 57565863 


54 


1390 


0 


2400 


29 . 54070381 


13 


2350 


1000 


4000 


21.93072188 


26 


1480 


0 


3900 


15 . 87080765 


20 


575 


0 


4200 


17.17841399 


3 


3200 


2500 


4000 


21 .00486097 


36 


760 


0 


2000 


47.1902545 


63 


3650 


0 


1100 


34.09124367 


44 


575 


1000 


4200 


56 .30447038 


11 


3300 


2500 


1200 


68 . 54997163 


46 


1565 


1000 


4000 


60.75108343 


48 


500 


0 


1000 


62.3335483 


51 


3450 


2500 


2100 


73.36103621 


35 


3740 


4500 


4500 


59 . 14895057 


23 


3600 


0 


500 


58 . 76346389 


10 


1500 


1000 


2500 


63.60363281 


43 


3300 


1000 


1000 


74.20451018 


18 


1450 


2500 


3800 


83.01091652 


27 


3300 


4500 


2500 


81.68791555 


49 


2480 


0 


1100 


45 . 51254274 


57 


2475 


2500 


4200 


74 .05848734 


39 


3550 


1000 


500 


89 .83687224 


42 


1620 


0 


1200 


61 . 06679397 


53 


2600 


1000 


2200 


66 .67038912 


1 


2400 


2500 


2200 


74 .05848734 


29 


2600 


0 


500 


83 .7322288 


2 


1600 


0 


500 


80 .96684485 


17 


2450 


1000 


1000 


71.82423312 


41 


2305 


4500 


4000 


89 . 08588313 


50 


1500 


4500 


4000 


111 . 141278 


58 


1550 


1000 


1100 


87.60219509 


33 


2640 


2500 


1100 


98.29833015 


24 


550 


2500 


4200 


87.38521723 


64 


625 


1000 


2200 


85.37050872 


47 


3500 


4500 


1200 


98.81618467 


55 


3475 


2500 


500 


107 . 5665614 


19 


3400 


4500 


500 


110 . 1665762 


5 


2500 


1000 


500 


96.09850519 


56 


500 


4500 


4000 


114 . 7626146 


32 


400 


4500 


2300 


132.4709759 


25 


2625 


4500 


2200 


99.42639298 


8 


640 


2500 


2400 


104.3943316 


12 


800 


0 


500 


95.50658187 


28 


660 


1000 


1000 


106.7833145 


62 


1500 


2500 


2400 


99.99999698 
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9 

34 

22 

14 

60 

52 

45 

30 

16 

6 

4 

61 

38 



2500 


4500 


500 


115 . 6789657 


1475 


1000 


500 


103 . 4813814 


1525 


2500 


1100 


106.5511138 


1575 


2500 


500 


119.3500945 


540 


1000 


500 


121.7141701 


700 


2500 


1000 


132.6538113 


2375 


2500 


500 


120 . 7684663 


1585 


4500 


2300 


116.6361357 


700 


4500 


1100 


139.5187601 


1400 


4500 


1100 


124.5815073 


750 


2500 


500 


142.2335988 


2500 


4500 


1000 


130 . 9359445 


1480 


4500 


500 


142.6119009 
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APPENDIX C. REGRESSION ANALYSIS OF COGNITIVE FACTORS 

AFFECTING MOVEMENT 



A. SURFACE PLOTS 

Surface plots of all possible combinations of independent variables were done in 
order to determine if there were any particular relationships between variables with re- 
spect to TMP which needed to be accounted for in the regression analysis. These are 
used in a fashion similar to scatter plots in simple regression. Note the linear or planar 
relationship between most combinations shown. 
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3-D PLOT OF TMP VERSUS RANGE AND SPEED 
CONSTANT LINE OF SIGHT = 0 










Figure 3. Surface Plot of TMP Versus SPD and RNG for No LOS 
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TMP IN INTERVAL SCALE 



3-D PLOT OF TMP VERSUS RANGE AND SPEED 
CONSTANT LINE OF SIGHT = 2 OR 3 




Figure 4. Surface Plot of TMP Versus SPD and RNG for 2-3 LOS 
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3-D PLOT OF TMP VERSUS RANGE AND SPEED 
CONSTANT LINE OF SIGHT = 4 OR MORE 




Figure 5. Surface Plot of TMP Versus SPD and RNG for 4- LOS 
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TMP IN INTERVAL SCALE 



3-D PLOT OF TMP VERSUS RANGE AND SPEED 
ALL LINES OF SIGHT 




Figure 6. Surface Plot of TMP Versus SPD and RNG for All LOS 
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TMP VERSUS SPEED AND NUMBER OF LINES OF SIGHT 

ALL RANGES 




Figure 7. Surface Plot of TMP Versus SPD and LOS for All RNG 
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TMP IN INTERVAL SCALE 



Figure 8 



TMP VERSUS RANGE AND NUMBER OF LINES OF SIGHT 

ALL RANGES 




Surface Plot of TMP Versus RNG and LOS for All SPD 
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B. REGRESSION USING GRAFSTAT 

MULTIPLE REGRESSION USING BACKWARD ELIMINATION 

TABLE OF COEFFICIENTS 

63 OBSERVATIONS R-SQUARED = 0.9048 STANDARD ERROR = 12.486 
5 VARIABLES ADJ R-SQUARED = 0.8983 



COEF 

INTERCEPT 
RNG3 
L0S3 
SPD 3 

(LOS3-SPD3 ) 



ESTIMATE 
1 . 1824F2 
-1.4744E-2 
1.1545E-2 
-1 . 9871E-2 
1.9839E-6 



STD ERR 
5 . 2563F0 
1.4786E-3 
1 .6888E-3 
1 . 7745F-3 
6 . 8170E-7 



T ST AT 
22.495 
-9.9715 
6.8366 
-11 . 198 
2 . 9102 



SIC LEVEL 
1.6653^-16 
3 . 4722E-14 
5 . 5107F-9 
6. 3838E-16 
5.1154E-3 



The assumptions of normality in the residuals and homogeneity of variance in the 
residuals are supported. The model is valid based on the standard error of the estimate 
and the extremely high R 2 . 



ANOVA RESULTING FROM THE REGRESSION 



ANALYSIS OF NORMAL DISTRIBUTION FIT 



DATA 

SELECTION 
X AXIS LABEL 
SAMPLE SIZE 
CENSORING 
FREQUENCIES 
EST . METHOD 
CONF METHOD 



RESIDUAL 

ALL 

RESIDUAL 

63 

NONE 

1 

MAXIMUM LIKELIHOOD 

ASYMPTOTIC NORMAL APPROXIMATION 



PARAMETER 

MU 

SIGMA 



ESTIMATE 

2 . 0510£”l3 
1 . 1980E1 



CONF. INTERVALS 
(95 PERCENT ) 
LOWER UPPER 
”2.9589 2.9589 

9.8878 14.072 



COVARIANCE MATRIX OF 
PARAMETER ESTIMATES 
MU SIGMA 

2.2781 0 

0 1.1391 





SAMPLE 


FITTED 


GOODNESS OF 


FIT 


MEAN : 


1.933 IE 13 


2 . 0510E”13 






STD DEV : 


1 . 2 0 7 6E1 


1 . 1980£1 


CHI-SQUARE 


2 .9566 


SKEWNESS i 


-3 . 1684E-2 


0.0000F0 


DEG FREED 


3 


KURTOSIS : 


2 . 6120E0 


3 . 0000E0 


SIGNIF 


0 .39836 








KOLM-SMIRN 


0 .060943 


PERCENTILES 


SAMPLE 


FITTED 


SIGNIF 


0 .97341 


5 : 


-20.102 


-1 . 9710E1 


CRAMER -V M 


0 .037255 


10 : 


-14.305 


-1 . 5355E1 


SIGNIF 


> .15 


25: 


-8.0305 


-8 . 0768E0 


ANDER-DARL 


0 .22113 


50: 


-1 . 0027 


1 . 2101E-6 


SIGNIF 


> .15 


75: 


10.122 


8 . 0768E0 






90 : 


14.321 


1 . 5355E1 







S2 



DENSITY FUNCTION 



CUMULATIVE DISTRIBUTION FUNCTION 







Figure 9. Residual Analysis of the Four Variable Regression 
9 5 : 19.153 1.9710E1 
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GOODNESS OF FIT ON THE RESIDUALS FROM THE REGRESSION 



CHI-SQUARE GOODNESS OF FIT TABLE 



LOWER 


UPPER 


OBS 




-INF. 


15.639 


6 


6 


15.639 


7 . 8194 


11 


10 


“ 7.8194 


0 


16 


15 


0 


7.8194 


11 


15 


7.8194 


15.639 


14 


10 


15.639 


+INF. 


5 


6 


TOTAL 




63 


63 



EXP 
. 0403 
. 149 
. 311 
.311 
.149 
.0403 

2.9566 



0 

0 

0 

'4 

3 

'1 



O-E 

, 040285 

,85093 

,68937 

,3107 

.8509 

, 0403 



((0-E)*2)*E 



0 

0 

0 

1 

1 

0 



00026868 

071344 

031039 

2136 

4612 

17916 
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DENSITY FUNCTION 



CUMULATIVE DISTRIBUTION FUNCTION 







Figure 10. Residual Analysis of the Five Variable Regression 
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C. REGRESSION USING SAS 



1. THE FIRST SAS MULTIPLE REGRESSION 



OPTIONS LINESIZE=80 ; 

DATA ONE ; „ 

CMS FILEDEF DD 1 DISK GDATA22 DATA A ; 



INFILE DD 1; 

INPUT RANG LOSS SPED TAMP ; 

PROC STEPWISE DATA=ONE; 

MODEL TAMP = PAJVC LOSS SPSS / BACfWAPL MAXR ; 



SAS 9:54 TUESDAY, JUNE 13, 1989 

BACKWARD ELIMINATION PROCEDURE FOR DEPENDENT VARIABLE TAM. 



STEP 0 


ALL VARIABLES ENTERED 


R SQUARE = 0. 
C(P) = H. 


89099524 

00000000 








DF 


SUM OF SQUARES 


MEAN SQUARE 


F 


PR0B>F 


REGRESSION 

ERROR 

TOTAL 


3 

59 

62 


84699.83472023 

10362.21574916 

95062.05046940 


28233 .27824008 
175.63077541 


160.75 


0.0001 




B 


VALUE 


STD ERROR 


TYPE II SS 


F 


PR0B>F 


INTCPT 

RANG 

LOSS 

SPED 


110 

-0 

0 

-0 


. 47735149 

! 01482470 0.00156910 
.01563357 0.00099485 
.01595209 0.00122666 


15677.32260814 

43371.36043090 

29702.31143979 


89.26 
246 .95 
169.12 


0 . 0001; 
0 .0001, 
0 . oooi 


BOUNDS 


ON 


CONDITION NUMBER: 


1.002224, 9.013743 





ALL VARIABLES IN THE MODEL ARE SIGNIFICANT AT THE 0.1000 LEVEL . 
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2. THE SECOND SAS MULTIPLE REGRESSION 



OPTIONS LINESIZE- 80; 

DATA ONE; 

CMS FILEDEF DD1 DISK CDATA22 DATA A; 

INFILE DD1; 

INPUT RANG LOSS SPED TAMP; 

RS - RANG* SPED; 

LS = L0SS*SPED; 

RL = RANG* LOSS; 

PROC STEPWISE DATA=0NE; 

MODEL TAMP - RANG LOSS SPED RS LS RL / BACKWARD MAXR ; 



SAS 9:55 TUESDAY, JUNE 13, 1989 

BACKWARD ELIMINATION PROCEDURE FOR DEPENDENT VARIABLE TAMP 



STEP 0 


ALL VARIABLES 


ENTERED 


R SQUARE = 0. 


90984757 












C(P) = 7. 


00000000 








DF SUM OF SQUARES 


MEAN SQUARE 


F 


PR0B>F 


REGRESSION 


6 86491 


. 97593599 


14415.32932267 


94 . 20 


0 . 0001 


ERROR 




56 8570. 


07453341 


153.03704524 






TOTAL 




62 95062 


.05046940 










B 


VALUE 


STD ERROR 


TYPE II SS 


F 


PR0B>F 


INTCPT 


109 


. 07985496 










RANG 


-0 


. 01017627 


0 . 00299509 


1766 . 66050270 


11.54 


0 .0013 


LOSS 


0 


.01337164 


0 .00249723 


4387 . 82533431 


28.67 


0 .0001 


SPED 


-0 


.01728175 


0 .00270126 


6263 . 83410262 


40 .93 


0 .0001 


RS 


-0 


.00000133 


0 .00000106 


240 . 98633319 


1 . 57 


0 .2147 


LS 


0 


. 00000203 


0 .00000068 


1373 . 72045367 


8 .98 


0 .0041 


RL 


-0 


.00000093 


0 . 00000086 


179.42161660 


1.17 


0 .2835 


BOUNDS 


ON 


CONDITION NUMBER : 7 


.239146, 206. 


9791 
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STEP 1 


VARIABLE RL 


REMOVED 


R SQUARE = 0 
C(P) = 6 


.90796016 

.17240643 






DF SUM OF SQUARES 


MEAN SQUARE 


F 


PROB>F 


REGRESSION 5 86312.55431939 

ERROR 57 8749.49615001 

TOTAL 62 95062.05046940 


17262 . 51086388 
153.49993246 


112.46 


0 .0001 




B VALUE 


STD ERROR 


TYPE II SS 


F 


PROB>F 


INTCPT 

RANG 

LOSS 

SPED 

RS 

LS 


112.32206019 
-0.01170732 
0.01136915 
-0.01705147 
-0 . 00000146 
0.00000206 


0 .00264429 
0.00168059 
0 .00269694 
0 . 00000105 
0 .00000068 


3008 . 87281840 
7024.89975456 
6136 . 04920158 
292.39366190 
1411 . 56803480 


19.60 
45 . 76 
39 . 97 
1.90 
9 .20 


0.0001 
0 .0001 
0 .0001 
0.1729 
0 . 0036 


BOUNDS 


ON CONDITION 


NUMBER : 


6.576496, 116 


.9728 




STEP 2 


VARIABLE RS 


REMOVED 


R SQUARE = 0 
C(P) = 6 


.90488434 
. 08301351 






DF SUM OF SQUARES 


MEAN SQUARE 


F 


PROB>F 


REGRESSION 4 86020.16065748 
ERROR 58 9041.88981191 
TOTAL 62 95062.05046940 


21505 . 04016437 
155.89465193 


137 .95 


0 . 0001 




B VALUE 


STD ERROR 


TYPE II SS 


F 


PROB>F 


INTCPT 

RANG 

LOSS 

SPED 

LS 


118 .23954426 
-0.01474358 
0 .01154539 
-0.01987087 
0.00000198 


0 . 00147857 
0 . 00168875 
0 . 00177449 
0 . 00000068 


15500 .71093605 
7286.43553484 
19548 . 66357433 
1320.32593725 


99.43 
46 .74 
125.40 
8.47 


0 . 0001 
0.0001 
0.0001 
0 . 0051 


BOUNDS 


ON CONDITION 


NUMBER : 


4.72315, 


45.35 





all variables in the model are significant AT THE 0.1000 LEVEL. 

SUMMARY OF BACKWARD ELIMINATION PROCEDURE FOR DEPENDENT VARIABLE 



STEP 

1 

2 



VARIABLE 

RMVED 

RL 

RS 



NUMBER PARTIAL 



IN 

5 

4 



R** 2 

0.0019 

0.0031 



MODEL 

R** 2 

0 . 9080 
0.9049 



C(P) 

6 . 17241 
6 . 08301 



1 .1724 
1.9048 



PROB>F 

0 . 2835 
0 .1729 
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1989 



SAS 9:55 TUESDAY, JUNE 13, 

MAXIMUM R- SQUARE IMPROVEMENT FOR DEPENDENT VARIABLE TAMP 



STEPS 


1 THROUGH 3 OMITTED 








STEP 4 


VARIABLE LS ENTERED 


R SQUARE 
C(P ) = 


= 0.90488434 
6 .08301351 




DF SUM OF SQUARES 


MEAN SQUARE 


F 


PROB>F 


REGRESSION 4 86020.16065748 
ERROR 58 9041.88981191 
TOTAL 62 95062.05046940 


21505 . 04016437 
155 . 89465193 


137.95 


0.0001 




B VALUE STD ERROR 


TYPE II SS 


F 


PROB>F 


INTCPT 

RANG 

LOSS 

SPED 

LS 


118 .23954426 

-0.01474358 0.00147857 

0.01154539 0.00168875 

-0.01987087 0.00177449 

0.00000198 0.00000068 


15500.71093605 

7286.43553484 

19548.66357433 

1320.32593725 


99.43 
46 .74 
125.40 
8.47 


0 . 0001 
0 . 0001 
0 . 0001 
0 .0051 


BOUNDS 


ON CONDITION NUMBER: 4 


.72315, 4 


5.35 





THE ABOVE MODEL IS THE BEST 4 VARIABLE MODEL FOUND. 



STEP 5 VARIABLE 





DF 


REGRESSION 


5 


ERROR 


57 


TOTAL 


62 



B VALUE 



RS ENTERED 



SUM OF SQUARES 

86312.55431939 
8749 . 49615001 
95062.05046940 

STD ERROR 



R SQUARE 
C(P) = 

MEAN SQUARE 

17262.51086388 

153.49993246 



TYPE II SS 



3008.87281840 
7024 .89975456 
6136.04920158 
292.39366190 
1411.56803480 

576496, 116 



= 0.90796016 
6 .17240643 

F PROB>F 

112.46 0.0001 



F PROB>F 



19.60 0.0001 

45.76 0.0001 

39.97 0.0001 

1.90 0.1729 

9.20 0.0036 

. 9728 



INTCPT 112.32206019 



RANG 

LOSS 

SPED 

RS 

LS 



-0 , 
0 , 
-0 , 
- 0 , 
0 , 



01170732 

01136915 

01705147 

00000146 

00000206 



0 

0 

0 

0 

0 



00264429 

00168059 

00269694 

00000105 

00000068 



BOUNDS ON CONDITION NUMBER: 



THE ABOVE MODEL IS THE BEST 5 VARIABLE MODEL FOUND . 
STEP 6 OMITTED 
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APPENDIX D. 



COMBINATORIAL OPTIMIZATION PROGRAM USING 
LAG RAN GIAN RELAXATION 



PROGRAM NTWRKS 

' V 'V ‘V /V Vv Vc Vr Vf Vr *V Vr V? * V V: Vc -,V 'V *V Vr Vr ‘V Vr Vr Vr *V “V iV iV Vr Vr ■?> Vr “V "sV ">V Vc Vr Vr Vc Vr Vc Vc /V * V ‘V * Y ?Y Vc * V * * V ‘V Vr Vr Vc "jV * V * V Vr Vr Vc Vc Vc Vc ' 



* 

* 

Vr 



* 

* 

* 

* 

* 

“V 



PROBLEM: 



>v* 

* 

* 

* 

* 

* 

* 

* 

* 

* 



AUTHORS: 



* SUBROUTINE: 



* VARIABLES: 


UVCD - 


sV 




*Y 


EDGES ■ 


V? 


XI, X2 ■ 


-,v 


OPTX - 


V? 


NODES 


V: 


S, T - 


* 


U, V - 


V- 


EPS - ] 


* 


TBAR - 


* 


PRED - 


y? 


D - D(’ 


* 


DMIN - 


v? 


EP - El 


* 


OUTDEG 


* 


ADJ - , 


“V 


LENGTH 


* 


TACDIF 


* 


TDMAX 


* 


OBJ - 1 




OBJMAX 


* 


LAM - : 




LAMM IN 


“V 


LAMM ID 


* 


LAMMAX 


* 


ENDLAM 


* 


TDLAMT 


iV 


SLOPE 



A TACTICAL UNIT MUST MOVE FROM START POINT A, THROUGH 
AN INTERMEDIATE POINT B, TO OBJECTIVE C IN NO MORE THAN 
TBAR TIME UNITS. EACH ARC/EDGE (U,V) IN THE NETWORK 
HAS A TRAVERSAL TIME TC(E) AND EXPOSURE/TACTICAL DIFFI- 
CULTY MOVEMENT ANTI -POTENTIAL OF TD(E). WE WISH TO FIND 
A PATH FROM A THROUGH B TO C WHICH HAS AN ACCEPTABLE 
TRAVERSAL TIME AND TACTICAL DIFFICULTY OR TO MINIMIZE 
TD SUBJECT TO TBAR. IT IS ASSUMED THAT THERE EXISTS AT 
LEAST 1 FEASIBLE PATH FROM A THROUGH B TO C. 

RICHARD MILLER, MAJOR, USMC 
CHARLES H. SHAW, CAPTAIN, USA 
IAN KEITH, CAPTAIN, USMC 
DYKSTR - DIJKSTRA' S ALGORITHM 

INPUT ARRAY OF EDGES (U,V),TIME COST TC(E), AND 
TACTICAL DIFFICULTY TD(E) . 

- NUMBER OF ARCS/EDGES (lOOOO MAXIMUM). 

- BOOLEAN VECTORS FROM X IN SUBROUTINE DIJKSTRA. 



NUMBER OF NODES, I.E. , VERTICES (1100 MAXIMUM). * 
START/STOP VERTICES FOR DIJKSTRA’ S ALGORITHM. * 

ALL ARBITRARY VERTICES. 

EPSILON, A SMALL NUMBER OR DIFFERENTIAL. * 

MAXIMUM TIME ALLOWED IN NETWORK (INPUT). * 

PREDECESSOR OF A NODE. 

) ’LABEL’ ARRAY FOR DIJKSTRA’ S ALGORITHM. * 

DMIN FOR DIJKSTRA’ S ALGORITHM. * 

TRY POINT ARRAY. 

- OUTDEGREE OF NODE TO COMPUTE ENTRY POINT ARRAY. * 
ADJACENT VERTICES ARRAY. 



EXPOSURE/TACTICAL DIFFICULTY ARRAY. 
LARGEST VALUE OF TACDIF( ). 

OBJECTIVE FUNCTION VALUE OF PI IN DIJKSTRA. 

OBJECTIVE FUNCTION MAXIMUM FOR P3. 
LAGRANGIAN MULTIPLIER IN SUBROUTINE DIJKSTRA. 
TEMPORARY MINIMUM OF LAM. 

MIDPOINT OF LAMM IN AND LAMMAX. 

TEMPORARY MAXIMUM OF LAM. 

BEST VALUE OF LAM. 

(TD + LAMBDA*T) ARRAY. 

SLOPE (TC*X - TBAR) IN SUBROUTINE DIJKSTRA 



7C 

* 

* 

Vc 

* 
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INTEGER UVCD( 10000,4) ,X1( 10000) ,X2( 10000) ,OPTX( 10000) 
INTEGER EP( 10001) ,OUTDEG( 10001) 

INTEGER ADJ( 10001) , LENGTH! 10001) ,TACDIF( 10001) 

INTEGER EDGES,I,J,K, NODES, A, B,C, TDMAX, L, TIME, U,V 

REAL TBAR , LAMM IN , LAMMID , LAMMAX , OB JMAX , EPS , P 1 , END LAM 
REAL SLOPE 1 , SLOPE2 , SLOPE3 , OB J1 , OB J2 , OB J3 

COMMON UVCD , EP , AD J , LENGTH , TACDIF , EDGES , NODES , TB AR , A , B , C 

PARAMETER (EPS = 0.0001) 

* SET FILEDEFS 

CALL EXCMS( ’FILEDEF 01 DISK PROJNET DATA B') 

CALL EXCMSC ’FILEDEF 02 DISK NETWORKS OUTPUT A') 

* INITIALIZE ARRAYS TO ZERO 

DO 5 I = 1, 10001 
OUTDEG(I) = 0 
EP( I ) = 0 
ADJ(I) = 0 
LENGTH! I) = 0 
TACDIF! I) = 0 
5 CONTINUE 

* READ IN THE DATA ONE ROW AT A TIME AND FIND 

* THE NUMBER OF NODES AND MAXIMUM VALUE OF TD(E) 

READ(1,*) NODES, A, B,C, TBAR 
I = 1 
TDMAX = 0 

10 READ! 1,*,END=20) (UVCD(I,J), J=l,4) 

TDMAX = MAX(UVCD( 1,4), TDMAX) 

1 = 1 + 1 
GO TO 10 

20 IF (NODES .GT. 1100) THEN 

STOP '*** ERROR, NUMBER OF NODES EXCEEDS 1000 ***' 
END IF 

* COMPUTE NUMBER OF EDGES IN THE NETWORK 

EDGES =1-1 

IF (EDGES .GT. 10000) THEN 

STOP '*** NUMBER OF EDGES EXCEEDS SIZE DEFINED ***' 
END IF 

* COMPUTE OUTDEGREE FOR EACH NODE 
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DO 30 I = 1, EDGES 

OUTDEG( UVCD( 1,1)) = OUTDEG( UVCD( 1,1)) + 1 
30 CONTINUE 

INITIALIZE ENTRY POINT ARRAY EP( ) 

EP( 1) = 1 

DO 40 I = 1, NODES 

EP( 1+1) = OUTDEG(I) 

EP( 1+1) = EP( 1+1) + EP( I ) 

40 CONTINUE 

COMPUTE THE HIERARCHICAL ADJACENCY LISTS FOR EACH NODE, EP( ) TO 
ADJ( ) AND LIST ALL THE NODE'S ADJACENT VERTICES AND CORRESPONDING 
LENGTH OR TIME COST AND EXPOSURE/TACTICAL DIFFICULTY 

K = 1 

DO 50 I = 1, NODES 

DO 51 J = 1, EDGES 

IF (UVCD( J, 1) .EQ. I) THEN 
ADJ(K) = UVCD( J,2) 

LENGTH(K) = UVCD(J,3) 

TACDIF(K) = UVCD(J,4) 

K = K + 1 
END IF 

IF (K . GT. 10000) THEN 

STOP '*** ADJACENCY LISTS EXCEEDS SIZE DEFINED ***' 
END IF 

51 CONTINUE 

50 CONTINUE 

PRINT *, 'INITIALIZATION COMPLETE AND BEGINNING SEARCH' 

FIND BOUNDS ON LAM 
LAMM IN = 0.0 

LAMMAX = (NODES - 1) * TDMAX 



L = 0 

CALL DYKSTR( LAMM IN , XI , SLOPE1 , OB J1 ) 

WRITE( 2 ,*) ' LAMMIN=' , LAMM IN, ' SLOPE 1= ’ ,SL0PE1 , ' 0BJ1=' ,0BJ1 
CALL DYKSTRC LAMMAX , OPTX , SLOPE 3 , OB J3 ) 

OBJMAX = OBJ3 

WRITE ( 2 ,*) ' LAMMAX= ' .LAMMAX, ' SLOPE3=’ ,SLOPE3, ' OBJ3=' ,OBJ3 
PRINT *,' INITIAL DIJKSTRA COMPLETE' 

70 LAMMID = ( LAMM IN + LAMMAX)/2.0 

IF ( ABS( LAMMAX - LAMMIN) . LE. EPS) GO TO 500 
L = L + 1 

PRINT * , ' ITERATION’ , L, ' STARTING* 
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CALL DYKSTRC LAMMID , X2 , SL0PE2 , OB J2 ) 

WRITE ( 2 ,*) 'LAMMID=' , LAMMID, ' SLOPE2=’ , SLOPE2 , ' OBJ2=' ,OBJ2 



IF (SLOPE2 . GT. 0) THEN 
LAMM IN = LAMMID 
END IF 

IF ( SL0PE2 . LE. 0) THEN 
LAMM AX = LAMMID 
IF ( 0BJ2 .GT. OBJMAX) THEN 
OBJMAX = OBJ2 
ENDLAM = LAMMID 
DO 80 I = 1, EDGES 

OPTX(I) = X2( I ) 

80 CONTINUE 

END IF 
END IF 
GO TO 70 

500 WRITEC2,*) 'WITH' ,L, ' ITERATIONS THE ' 

WRITE( 2,85) OBJMAX 

85 FORMATC MAXIMUM OBJECTIVE FUNCTION VALUE FOR P2 IS ' ,F10. 3) 
TIME = 0 
PI = 0. 0 

DO 90 I = 1, EDGES 

TIME = OPTX(I) * UVCD( 1,3) + TIME 
PI = OPTX(I) * UVCD( 1,4) + PI 
90 CONTINUE 

WRITE( 2 ,*) 'THE SOLUTION FOR (TD*X) IS ' ,P1 

WRITEC2,*) 'EPSILON IS ' , PI - OBJMAX,', TIME ON THE PATH IS ’ 
WRITE( 2 ,*) ’AND THE LAST VALUE OF LAM IS ’, ENDLAM 
WRITE(2,*) ’THE BEST PATH FROM NODE * , A , ’ , THROUGH NODE ’ ,B , ' 
WRITEC 2 ,*) ’TO OBJECTIVE NODE ’ ,C, ' IS ALONG PATH: ' 

WRITE( 2,95) A 
95 FORMAT( 'NODE ',15) 

U = A 
550 I = 1 

600 IF ((U . EQ. UVCD( 1,1)) .AND. (OPTX(I) . EQ. 1)) THEN 
WRITEC 2, 95) UVCDC 1,2) 

IF (C .EQ. UVCDC I, 2)) STOP 'THE END' 

U = UVCDC I, 2) 

GO TO 550 

ELSE 

1 = 1 + 1 
GO TO 600 
END IF 

END 
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SUBROUTINE DYKSTR( LAM , X , SLOPE , OB J) 



INTEGER PRED( 1100) ,1 , J,K,S,T,U,V, VMIN,X( 10000) 

INTEGER UVCD( 10000,4) ,EP( 10001) ,ADJ( 10001) , LENGTH ( 10001) 
INTEGER TACDIF( 10001) , EDGES , NODES ,A,B ,C 

REAL D( 1100) , DMIN, SLOPE, OBJ, LAM, TDLAMT( 10001) , TINF,TBAR 

COMMON UVCD , EP , AD J , LENGTH , TACD IF , EDGES , NODES , TBAR , A , B , C 

PARAMETER (TINF = 1.0E15) 

SET TDLAMTCO FOR NEW LAM 

DO 300 I = 1, EDGES + 1 

TDLAMT(I) = TACDIF(I) + LAM * LENGTH(I) 

300 CONTINUE 

RESET OPTIMAL PATH FOR NEXT SEGMENT 

DO 305 I = 1, EDGES 
X(I) = 0 
305 CONTINUE 

DIJKSTRA'S ALGORITHM FOR SHORTEST PATH FROM S TO ALL NODES 

S = A 
T = B 

INITIALIZE THE SHORTEST PATH AND PREDECESSOR ARRAYS 

400 DO 310 I = 1, NODES 
PRED(I) = 0 
D(I) = TINF 
310 CONTINUE 

INITIALIZE START VERTEX, S 

D(S) = 0. 0 
PRED(S) = S 

DO 320 I = 1, NODES 
DMIN = TINF 
DO 330 J = 1, NODES 

IF ((PRED(J).GE. 0) .AND. (D( J) . LT. DMIN) ) THEN 
DMIN = D(J) 

VMIN = J 
END IF 

330 CONTINUE 

PRED(VMIN) = -PRED(VMIN) 
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DO 340 K = EP(VMIN) , EP(VMIN+1)-1 

IF ((D(VMIN) + TDLAMT(K)) . LT. D( ADJ(K) ) ) THEN 
D(ADJ(K)) = D(VMIN) + TDLAMT(K) 

PRED(ADJ(K) ) = VMIN 
END IF 

340 CONTINUE 

320 CONTINUE 

* CONVERT PRED() TO PATH X( ) 

V = T 
360 I = 1 

370 IF ((-PRED(V) . EQ. UVCD(I,1)) .AND. (V . EQ. UVCD(I,2))) THEN 
X(I) = 1 
V = -PRED(V) 

GO TO 360 

ELSE 

IF (V .EQ. -PRED(S)) GO TO 380 
1 = 1 + 1 
GO TO 370 
END IF 

380 IF (S .EQ. B) GO TO 390 

* NOW GET SHORTEST PATH FROM B TO C 

S = B 
T = C 
GO TO 400 

* COMPUTE OBJECTIVE FUNCTION VALUE AND SLOPE FOR LAM 

390 OBJ = 0. 0 
SLOPE = 0. 0 
DO 450 I = 1, EDGES 

OBJ = (UVCD( 1,4) + LAM * UVCD(I,3)) * X(I) + OBJ 
SLOPE = UVCD( 1,3) * X(I) + SLOPE 
450 CONTINUE 

OBJ = OBJ - LAM * TBAR 
SLOPE = SLOPE - TBAR 

END 
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APPENDIX E. EXAMPLE MODEL OUTPUT 



The purpose of these example output plots is to demonstrate proper functioning of 
the model in terms of tactical movement decision logic and some sensitivity analysis of 
the cognitive time scale values. Plots are appropriately labelled with all normal Tactical 
Optimum Path (TOP) output plots being expected value model runs unless otherwise 
noted. High risk averse plots use an input value of 30 minutes and extreme risk averse 
plots use an input value of 60 minutes rather than the expected value of 5.5 minutes. 

These verification model runs are based on the Lauterbach map sheet near Fulda in 
the central region of The Federal Republic of Germany. The base map is one version 
of the terrain relief map available in the CAM MS. The NATO vehicle represented is 
an M1A1 Main Battle Tank moving as part of a deliberate attack against the specified 
Warsaw Pact systems shown in a deliberate defense. 

The sequence of numbers shown on the mapping represents the path points in the 
order specified by the user. The continuous line connecting these points is the resulting 
optimum path from the optimization. The other irregular polygon shapes represent the 
user input area of operations which appropriately restricts the network searched during 
the optimization. The set of optimum path points and large amounts of audit data are 
available in data files for use directly in combat models and analysis rather than output 
to a screen or map. 
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Figure 13. CAMMS Output Mapping For Set 4 
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Figure 18. TOP Output Mapping For Set 6 
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Figure 19. CAMMS Output Mapping For Set 15 
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